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ABSTRACT 

Complementing the observations on the eastern part of galaxy NGC 6946 presented in a previous work, we 
report high spatial resolution observations of Giant Molecular Clouds (GMCs) in the nearby spiral galax¬ 
ies M101 and NGC 628 obtained with the Combined Array for Research in Millimeter-wave Astronomy 
(CARMA). We observed CO(l —J 0) over regions with active star formation extending from 2 kpc to 15 kpc 
galactocentric radius. Higher resolution observations of CO(2 —J 1) toward the brightest regions observed in 
CO(l —> 0) have allowed us to resolve some of the largest GMCs. We have recovered short-spacing u-v com¬ 
ponents by using single dish observations from the Nobeyama 45m and IRAM 30m telescopes. Using the 
automated CPROPS algorithm we identified 112 CO complexes in the COG —J 0) maps and 144 GMCs in the 
CO(2 —> 1) maps. Using a Bayesian fitting approach, we generate scaling relations for the sizes, line widths, 
and virial masses of the structures identified in this work. We do not find evidence for a tight power law rela¬ 
tion between size and line width, although the limited dynamic range in cloud size remains a clear issue in our 
analysis. Additionally, we use a Bayesian approach to analyze the scaling relation between the star formation 
and molecular gas surface density, known as the Kennicutt-Schmidt relation. When we perform our analysis 
using the boundaries of the structures identified by CPROPS, we find that the distribution of slopes are broadly 
distributed, mainly due to the limited dynamic range of our measured E H 2 - In the case of the CO(l —J 0) 
complexes, the slope distributions are most consistent with super-linear relations, although sub-linear relations 
cannot be excluded for NGC 628 and NGC 6946. The GMCs from higher resolution CO(2 —> 1) maps follow 
a similar behavior, but with larger scatter. As a complementary study, we use the Bayesian approach to analyze 
the Kennicutt-Schmidt relation for a uniform grid covering the areas surveyed, and with Xh 2 non-detections 
included in the analysis. The distributions of slopes is consistent with sub-linear relations for NGC 6946 and 
NGC 628, but is less constrained for M101. This picture is preserved after a 24/rm background component is 
subtracted from the X$fr measurements. On-arm regions tend to have higher star formation rates than inter-arm 
regions. Similar to what we find in our study of the eastern part of NGC 6946, in M101 we find regions where 
the star formation efficiency (SFE) shows marked peaks at specific galoctocentric radii. On the other hand, the 
distribution of SFE in NGC 628 is more contiguous. We hypothesize that differences in the distribution of SFE 
may be indicative of different processes driving the spiral structure. 

Subject headings: galaxies: ISM — stars: formation — ISM: molecules — galaxies: individual (NGC 6946, 

NGC 628, M101) 


1. INTRODUCTION 

The formation of stars is the final stage of many physical 
processes acting on the complex structure of the ISM at dif¬ 
ferent temporal and spatial scales. As a consequence, a com¬ 
plete understanding of the processes responsible for how stars 
form out of the ISM remains a challenging task. However, it is 
observed that most of the massive star formation in the Milky 
Way (MW) is associated with GMCs, and the corrrelation be¬ 
tween star formation and molecular mass has been observed 
to hold up to much larger sizes, even ~ 1 kpc (Leroy et al. 


|2013 Bigiel et al.|20 08 Bigie l et al.|201 l| i. Therefore, study 
ing GMCs in different environments seems to be a crucial tool 
to connect the processes that trigger and regulate the star for¬ 
mation at large scales with those playing an important role 
for star formation within the clouds, allowing us to establish 
a unified picture of the star formation process across the full 
range of sizes observed. 

Due to observational limitations, most of our understanding 


of the physical processes involved in the creation of molecu¬ 
lar clouds, their subsequent fragmentation into smaller and 
denser structures such as clumps and cores, and finally the 
birth of clusters and individual stars, has come from the de¬ 
tailed exploration of individual star-forming regions in the 
Milky Way. On the other hand, extragalactic studies have 
focused mainly on the integrated star formation over regions 
commonly larger than 100 pc, or even over the entire galactic 
disk. Resolving structures close to the scales where GMCs 
form (~100 pc) in nearby galaxies is still observationally 
challenging. Thus, high resolution CO surveys of galaxies 
(|Sakamoto et al.||1999| |Helfer et al.||2003)> have focused pri¬ 
marily on th e centr al regions. Single-dish observations ( Kuno| 
|et al. 2007[ [Leroy et al.||2009l > sometimes offer a complete 
mapping of the disk, but the resulting resolutions have been 


inadequate to resolve GMCs. A complete census of GMCs 
generally requires a dedicated observing program over the 
disk of a single galaxy, such as has been conducted for M33 
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(Engargiola et al. 

2003 

Schruba et al. 2010), the Magel- 

lanic Clouds (NAN 

TEN 

surveys, Pukui et al. 2008), and M51 


In the past few years, numerical simulations have suggested 
that the formation and evolution of molecular clouds may de¬ 
pend on the dynamical process that drives the spiral arm struc¬ 
ture ( Dobbs & Bonnell|2006| |Dobbs & Bonnell|2008[ |Wada| 


et al. 201 1[ >. In the case of grand-design spiral galaxies, mas¬ 

sive molecular clouds are thought to be formed by the passing 
of spiral density waves, inducing the formation of spurs down¬ 
stream of the spiral arms due to the interaction between the 
ISM and the s tellar spiral arm (Kim & Qstriker||2002 Wada 


& Koda|2004|). Thus, an offset between the spir al shock and 


the potential minimum is expected (e. g. M51, |Louie et al. 
2013} >. On the other hand, in galaxies with a multi-arm spi 


ral structure, simulations find that the gaseous spiral structure 
coincides with the potential minimum over several crossing 
times. In this case, the collision or merging between spiral 
arms leads to the formation of the most massive structures, in¬ 
creasing the star formation in those overdense regions. There¬ 
fore, high resolution observations of the molecular gas in spi¬ 
ral galaxies with different morphologies can help disentangle 
the real processes involved in the formation of dense gas, and 
the subsequent formation of stars. 

1.1. Molecular Cloud Scaling Relations 

Scaling relations between properties of resolved structures 
have been useful tools to decipher the physical conditions 
present in the interstellar medium (ISM). Studies conducted 
in the MW showed that resolved properties such as velocity 
dispersion, luminosity and size of GMCs were correlated, fol¬ 


lowing empirical relations known as Larson’s laws (Larson 
1981; Solomon et al .|1987j >. The scaling relation between the 


line width and the size, commonly referred as the first Lar¬ 
son’s law, has been used to investigate the character of the tur¬ 
bulence in the ISM. Nevertheless, the dominant source of tur¬ 
bulence in the ISM, in particular whether it is internally or ex¬ 


ternally driven, remains controversial (see McKee & Ostriker 
|2QQ7| >. The second and third Larson’s laws dictate that the 
GMC’s are roughly self-gravitating and have approximately 
similar surface densities respectively. Although Larson’s re¬ 
lations have been shown to hold across a wide range of extra- 


galactic environments (Blitz et al. 2007 Bolatto et al. 2008), 

recent galactic ( 

Heyer et al. 2009) and extragalactic studies 

(Wong et al. 20 

1, Hughes et al.|2013 Colombo et al. J2014) 


ally find a weak correlation between the size and line width of 
the identified structures, while inferring genuine differences 
in their physical properties such as velocity dispersion, CO 
peak brightness and mass surface densities. 

Although Larson’s laws originally relate properties such as 
size, line width and luminosity, an extra relation can be built 
assuming virialization of the clouds. Virial masses derived 
from the size and line width can be compared to the luminos¬ 
ity of the GMCs, thus defining a scaling factor that relates the 
intensity of the CO emission with the IL column density, the 
Xcx) factor. This conversion factor is thought to depend on the 
metallicity and the local radiation field acting on the regions 
probed by the CO observations ( Leroy et al.|201 1| |Sandstrom| 


|et al.||2013| |Bolatto et al.||2013 1 . The current generation of 
radio facilities has provided CO maps of galaxies at a suf¬ 
ficient spatial resolution to resolve individual GMCs. Thus, 
the Xco factor can be inferred using the virial mass method. 
The average values thus obtained for the conversion factor in 


nearby spiral galaxies are within a factor of two of the MW 
value (Donovan Meyer et ak 2012; Donovan Meyer et al. 
|2013| [ Rebolledo et al.l|2012| ). However, the Xco factor can 


be estimated using other methods. For example, by assum¬ 
ing a constant H 2 depletion time in dwarf galaxies, |Schruba| 
[et al.| ( ]2012) > find Xco factor values more than one order of 
magnitude higher than spiral galaxies with solar metallicity. 
By assuming a power-law relation between the star form ation 
rate and the molecular gas surface density (see Section [L2| ), 
Blanc et al. ( |2013| ) measure the dependence of the Acofac¬ 
tor with the galactocentric radius in NGC 628. They find that 
the Xco factor increases with radius, which added to the ob¬ 
served metallicity gradient across the disk, implies that lower 
metallicity regions coincides with larger values of Xco ■ Alter¬ 
natively, optically thin dust emissi on can be used as an inde- 
pendent tracer of the gas |Leroy et al. |2011 , Sandstrom et al. 


is estimated from 

<5 g td s dust- where 


2013J. The surface density of the gas E gas 
the dust surface density Edust through E gas 
Ltd is the gas to dust ratio. To constrain L td , E~ gas is indepen¬ 
dently obtained from the Eh 2 = Ac q/cq and Ehi maps. For a 
set of galaxies in the Local Group, Leroy et al. ( |201 1) > solved 
simultaneously for Xco and Lid over regions witfTsize scales 
~ 1 kpc 2 . For galaxies with higher metallicity, they found 
Aco ~ 1-4.5 x 10 20 cm _2 (K km s -1 ) -1 , similar to MW val¬ 
ues. On the other hand, for systems with lower metallicity, 
they found larger Aco by a factor of 10 compared to the con¬ 
version factor obtained for high metallicity regions, similar to 
the conclusion obtained by|Schruba et aI7| ( |2012[ >. Addition¬ 
ally, |Sandstromerakl((20T3fToun3^3epression in Aco in the 
central part of the disk of some galaxies in their sample. 

1.2. Star Formation Relations 

Another relation used to study the properties of the ISM is 
the empirical power-law relation that connects the star forma¬ 
tion surface density (Esfr) and the gas surface density (E gas ) 
available to form stars, i. e., the relation given by Esfr oc E gas 
( [Kennicutt & Evans||2012| and references therein). This re¬ 
lationship, known as the Kennicutt-Schmidt law (abbreviated 
as K-S henceforth), has been the subject of several works that 


Evans et al.||2009[ Heiderman et al.||2010 Gutermuth et al. 

2011) and nearby galaxies (Bigiel et al. 

2008; 

Liu et al. 2011 

Rahman et al.[2011, Leroy et al. 2013 

Momose et al. 2013). 


lar gas content and star formation rates across a wide dynamic 
range, the intrinsic form of the K-S law has remained elusive. 
The different slopes reported in previous studies mostly fall 
into two groups: those that report linear relations (a ~ 1), and 
those that find a super-linear relations (a > 1) in the K-S law. 
Several explanations have been proposed to explain the wide 
range of slopes (~ 1-3) found in extragalactic studies, e. g., 
biases in the tracer used to estimate the star formation rate 
or the molecular gas, biases in the fitting procedure used to 
find the power-law, variations in the fraction of dense gas, and 
the size scales over which averaging is performed (|Liu et alT 


|2011[ [Calzetti et ~aT7| [2012 ; Mo mose et al.||20l3] ). More re- 
cently, it has been proposed that diffuse emission in the maps 
used to estimate the star formation rates (FUV, 24pm and Ho) 
may be the reason for the divergent results. [Liu et al. (|201 lj), 
for example, were able to reconcile the linear relation found 


by Bigiel et al. 

2008) 

Kennicutt et al. 

(2007) 


in M51a by subtracting or preserving 
a diffuse background component in the star formation tracers. 
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FlO. 1.— CO(l —> 0) integrated intensity contours overlaid on 3.6 //m map from SINGS for the northern arms of NGC 628. Contours begin at 3.56 K km s 1 
and are spaced by 2x3.56 K km s -1 . The noise at the center of the map is 1.5 K km s 1 . The dashed black line illustrates the region where full gain correction 
was made. 


TABLE 1 

Properties of sample galaxies. 


TABLE 2 

CARMA DATA PROPERTIES. 


Property 

NGC 6946 

NGC 628 

M101 

Property 

NGC 6946 a 

NGC 628 

M101 

Morph. a 

SABcd 

SA(s)c 

SAB(rs)cd 


CO(l —¥ 0) maps 



R.A. (J2000) a 

20:34:52.3 

01:36:41.7 

14:03:12.5 

Sensitivity (K) 

0.40 

0.42 

0.37 

Deck (J2000 a 

60:09:14 

15:47:01 

54:20:56 

Velocity Resolution (km s _ 

*) 2.5 

2.5 

2.5 

Distance (Mpc) 

5.5 b 

7.3 d 

7.4 e 

Angular Resolution (") 

5.2 x 5.0 

4.5 x 4.4 

4.4 x 4.2 

Incl. (°) c 

33 

7 

18 

Linear Resolution (pc) b 

136 

157 

154 

PA. (°) c 

243 

20 

39 


CO(2 —»■ 1) maps 



Notes. 




Sensitivity (K) 

0.53 

0.31 

0.29 


a NAS A/IPAC Extragalactic Database (NED), 
b 


Velocity Resolution (km s 1 ) 


2.5 


2.5 


2.5 


Tully (1988 


Angular Resolution (") 

2.5 x 2.0 

2.0 x 1.8 

2.0 x 1.9 

Walter et al. 

12008 . 

Lineal* Resolution (pc) b 

60 

67 

70 


Sharina et al. 


Ferrarese et al. 


1 996 . 

fZooo} 


Although the correlation between gas and star formation 
rate (SFR) surface densities has been observed to be tight 
when averaged over substantial areas of galaxies, recent high 
resolution observations of molecular gas in nearby galaxies 
have r evealed a poorer cor relation bet ween these t wo quan¬ 
tities (Schru ba et al.|[20T0 |Rebolledo et al.||2012[ hereafter 
Paper I). These recent studies concluded that the evolution of 
individual regions is the main source of the increased scatter 
observed at small scales. Similarly, Calzetti et al.]( [2012[ ) find 
that the stochastic sampling of the molecular cloud mass func¬ 
tion determines the slope and the scatter of their simulated K- 
S relation. They find that at scales below ~ 1 kpc the K-S 
relation is super-linear, while at scales ~ 1-2 kpc the relation 
becomes linear as the cloud mass function is fully sampled. 

We have recently presented a study of the properties of the 
molecular gas in the eastern part of the nearby galaxy NGC 
6946 (Paper I). We found that the properties of the clouds 
follow relations similar to those found previously for extra- 
galactic clouds. In particular, the trends exhibited by the 
CO(l —>■ 0) complexes presented in that study were consis- 


Notes. 

a Properties of NGC 6946 as in Paper I. 
b Linear resolution is calculated from \/b u 


, x b m 


tent with those found in the center of NGC 6946 by Donovan 
Meyer et al. (2012 1 (hereafter DM12), despite the differences 


in resolution and cloud identification algorithms; the excep¬ 
tion is the set of GMCs located within 400 pc of the galactic 
center with large velocity dispersions (> 10 km s -1 ). Addi¬ 
tionally, our virial mass-CO luminosity relation is consistent 
with the relation found by DM12 for the central part of NGC 
6946 and is consistent with our choice of the CO to H 2 con¬ 
version factor Xco=2 x 10 20 cm _2 (K km s -1 ) _I . 

In this paper, we compare the properties of molecular struc¬ 
tures found in NGC 6946 with those located in sub-regions of 
two other nearby spiral galaxies: M101 and NGC 628. We 
have continued our observational strategy previously used for 
NGC 6946 (Paper I) which consists of starting with large- 
scale, lower resolution CO(l -A 0) maps, and following up 
with higher resolution CO(2 —> 1) imaging on smaller areas. 
We select regions in the disks with active star formation, and 
near the transition radius from the Ho-dominated inner galac- 
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Fig. 2.— CO(l —> 0) integrated intensity contours overlaid on 3.6 //m map from SINGS for the regions observed in M101. Contours begin at 5.1 K km s 1 and 
are space d by 2x5.1 K km s . The noise at the center of the maps is 2.5 K km s 1 . The blue circles indicates the regions with enhanced star formation reported 
in Section l53l 


tic region to the Hi-dominated outer region. As with the 
galaxy NGC 6946, these two galaxies represent excellent tar¬ 
gets because of their proximity, high CO surface brightness, 
low inclination, and the availability of high-quality datasets 
at different wavelengths 


ter et al 
and multibam 


including HI THINGS maps (Wal- 
.|2008b, GALEX UV imaging ( |Gil de Paz et al.|2007 i 
ltiband Spitzer imaging from SINGS (Kennicutt et al. 


[2003) . 

We present our study as follows: in Section [2] we describe 
our observations of NGC 6946, M101 and NGC 628 using 
CARMA, and we describe the archival data at several wave¬ 
lengths that we include in our analysis. In Section[3]we sum¬ 
marize the technique used to identify GMCs and to measure 
their physical properties. In Section[4]we present some statis¬ 
tics of the measured properties, and we introduce the Bayesian 
inference method used to study the different ISM scaling rela¬ 
tions found in our sample of galaxies. In Section[5]we discuss 
the implications of our results, and we investigate whether the 
star formation properties of the clouds differ betwee n on- arm 
and inter-arm regions or between galaxies. Section [53] dis¬ 
cusses the limitations of our analysis. In Section [6] we sum¬ 
marize the work presented in this paper. 


2. DATA 

2.1. CARMA observations 
2.1.1. 3 mm 

We performed high spatial resolution observations of 
CO(l -> 0) for the selected regions in NGC 6946, NGC 628 
and M101 from July of 2009 to September 2012. We used the 
Combined Array for Research in Millimeter-wave Astronomy 
(CARMA) in E, D and C array configurations, which have 
baselines of 8.5-66 meters, 11-148 meters, and 30-350 me¬ 
ters respectively. In this section we describe the details of the 
NGC 628 and M101 maps, whereas we refer the reader to 
Paper I for the observation strategy and details of NGC 6946 


maps. TableHlshows the basic parameters of NGC 6946, NGC 
628 and MIOl. 


For NGC 628, we observed CO(l —I 0) in a 53-point mo¬ 
saic area, which covers 5.5x2.7 arcmin * 2 , corresponding to 
a physical scale of 11.7x5.7 kpc 2 at a distance of 7.3 Mpc 
(calculated from the luminosity of the three brightest blue su- 
Sharina et ah] 19961. 


pergiant stars, 
to have the C(5TT 


The correlator was set 
0) line in the upper side band (USB). 


We have placed three overlapping 62 MHz bands to cover 
the line. One band was centered in the CO(l —I 0) line fre¬ 
quency, while the other two were placed with an offset of 15 
MHz from the rest frequency, achieving a total velocity cov¬ 
erage of 240 km s _1 and a channel width of 0.4 km s -1 . A 
500 MHz wide-band was placed outside the spectral windows 
for calibration. We observed every pointing for 30 seconds, 
yielding a 26.5-minute observation time per cycle. At the be¬ 
ginning of the track we observed 3C454.3 and 3C84 as the 
passband calibrators, and at the end of every cycle, we ob¬ 
served 0108+015 as the gain calibrator. Calibration, imaging 
and deconvolution were performed using standard procedures 
of the MIRIAD software package. The clean CO(l -+ 0) map 
has a resolution of 4"50 x 4"37, a er rms of 0.423 K, a chan¬ 
nel size of 2.5 km s -1 , and a pixel size of 1". Figur e [T] sho ws 
the CO(l —> 0) map of the region observed (Section |2.3.3| ex- 
plains our approach to generate moment maps) overlaid on a 
3.6 //m image from SINGS (Kennicutt et al.|2003jl. 

In the case of M101, we observed CO(l -4 0) in two differ¬ 
ent areas of the southern region. The main area was covered 
with a 53-point mosaic, corresponding to a physical scale of 
13.7x5 kpc 2 at a distance of 7.4 Mpc (calculated from the lu¬ 
minosity of the Cepheids,[Ferrarese et a l.|2000| l. Additionally, 
we have observed a small region in the western region of the 


disk with a 7-point mosaic, corresponding to a physical scale 
of 1 kpc 2 . The correlator setup is the same as for the NGC 628 
observations. We have used 1153+495 as the gain calibrator. 
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Fig. 3.— CO(l —» 0) integrated intensity map of the region observed in NGC 628. Contours are as in FigurejTj Black circles indicate the regions where we 
performed CO(2 —> 1) observations. 



RA offset (arcsec) 

Fig. 4.— COll > 0) integrated intensity map of the regions observed in M101. Contours are as in Figure[2] Black circles indicate the regions where we 
performed CO(2 —> 1) observations. 


and 3C279 and 3C273 as passband calibrators. Every point¬ 
ing was observed for 40 seconds, but in this case the mosaic 
was separated into two cycles of 30 pointings, yielding a 20- 
minute observation time per cycle. The clean CO(l -4 0) map 
has a resolution of 4"37 x 4"16, a <7 rms of 0.370 K, a chan¬ 
nel size of 2.5 km s -1 , and a pixel size of 1". Figure [5] shows 
the CO( 1 —> 0) map of the regions observed overlaid on a on 
a 3.6 fi m image from SINGS. Table [2] summarizes the basic 
properties of our data cubes. 

2.1.2. 1 mm 

In order to resolve the largest GMCs in the galaxies, we 
performed high resolution observations in D array towards 
some CO complexes using the CO(2 —> 1) line. We selected 
10 regions in each galaxy based on the highest CO(l —> 0) in¬ 
tegrated intensity in the area covered by 3 mm observations. 
Figures[3]and|4]illustrate the regions targeted with CO(2 —> 1) 
observations in NGC 628 and M101 respectively. The ob¬ 


servations were taken between October and November 2011. 
For both galaxies, we set the correlator to have a 125 MHz 
band centered on the rest frequency of the CO(2 -+ 1) line 
in the LSB, and seven 500 MHz wide bands to observe the 
continuum at 1 mm. This correlator configuration yields a 
velocity coverage of 160 km s -1 and a channel width of 0.5 
km s -1 . Every pointing was observed for 90 seconds, which 
yields a total cycle time of 15 minutes. For NGC 628 we 
observed the gain calibrator 0108+015 after every cycle, and 
3C84 was used as passband calibrator, while for M101 we ob¬ 
served 1642+689 and 3C279 as gain and passband calibrators 
respectively. The clean CO(2 -+ 1) maps for NGC 628 have 
a 2.5 km s -1 velocity resolution, and an angular resolution of 
2"03x 1 "80 corresponding to ~ 70 pc. The rms noise of the 
clean maps is 0.312 K. For M101, the clean maps have a 2.5 
kms' 1 velocity resolution, and a resolution of 2"01xl"87 
corresponding to ~ 72 pc. The rms noise of the clean maps is 
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0.287 K. The pixel size of the maps is 0.5". The CO(2 —> 1) 
maps are illustrated in Figures[5]and[6]for NGC 628 and M101 
respectively. 

Through this paper, we convert the CO(2 —► 1) maps to 
CO(l —> 0) assuming a fixed /co(2 1) to /co(l —> 0) ra¬ 

tio of 1 across the regions observed, in order to be consistent 
with the work presented in Paper I. Although previous studies 
have used fixed values of this quantity (e. g., 0.8, Leroy et al. 


2008 henceforth L08), this ratio has been observed to vary 


from 0.6 to 1.0 in different parts of the disk of nearby spiral 
galaxies (e. g„ HERACLES, ILeroy et al.||2009| M 51,|Kodal 
jet al.|2012| . -- 


2.2. Single dish maps 
2.2.1. NRO 45m CO( 1 -t 0) map 

Extended flux is recovered by merging our CO(l —► 0) 
CARMA data with single dish maps from the Nobeyama 45- 
meter single dish telescope (NRO 45m). The galaxy NGC 
628 is part of the CARMA and NObeyama Nearby-galaxies 
(CANON) CO(l-O) Survey project, which combines observa¬ 
tions from the CARMA and NRO 45m telescopes to resolve 
GMCs in disks of nearby galaxies (Koda et al., in prepara¬ 
tion). The CO( 1 —>• 0) map of galaxy NGC 628 was obtained 
using the Beam Array Receiver System (BEARS) instrument. 
The FWHM of the 45 meter dish is 15" at the CO(l —» 0) rest 
frequency which is degraded to 19."7 after regridding. The 
rms noise of the CO(l —t 0) single dish map is 0.06 K (0.3 
Jy beam -1 ) in a channel width of 2.54 km s -1 , which corre¬ 
sponds to <t(Eh,) ~ 7 Mg pc -2 assuming a line width of 20 
km s -1 and Xco = 2 x 10 20 cm -2 (K km s -1 ) -1 . 

For the main mosaic observed in M101, we use NRO 45m 
observations from |Kuno et al.| ( j2007) > to recover the extended 
flux, and we refer the reader to that paper for the observation 
details. The rms noise in the maps is 0.07 K, and the FWHM 
of the 45 meter dish is 15" before regridding at the CO( 1 —► 0) 
rest frequency. The velocity resolution of the map is 5 km s -1 . 
On the other hand, since the area covered by the Kuno et al. 
(20071 observations does not overlap with the small region 
located in the western arm of Ml01 covered with a 7-pointing 
mosaic, we have used a CO(2 — » 1) single dish map from the 
HERACLES project (see|2.2.2|) to combine with our CARMA 
map instead, assuming a line ratio of 1 as noted in Section 
[ITU 


2.2.2. Heracles CO{ 2 —> 1) map 

Following the same approach we adopted in Paper I, we 
have combined our CO(2 —> 1) CARMA data with single dish 
maps from the HERACLES project (Leroy et al.|2009). HER¬ 
ACLES observed the CO(2 —> 1) line toward the full opti¬ 
cal disk of 18 nearby galaxies with the HERA receiver at 
the IRAM 30m telescope, producing maps with 13" and 2.6 
km s -1 of angular and velocity resolution respectively. The 
er rms in the data cube is ~ 20 mK, which yields an inte¬ 
grated intensity map with a noise level of ct(Eh,) ~ 2 M 0 pc -2 
assuming a velocity width of 20 km s -1 and Xco = 2 x 
10 20 cm -2 (Kkm s -1 ) -1 . 


2.3. CARMA-Single dish merging procedure 
2.3.1. Combination of CO( 1 —> 0) data 

We have used the approach introduced by Koda et al.|(j2011| 
to combine our CO(l —> 0) maps from CARMA with the 
single dish maps from NRO 45m single dish telescope. A 


detailed description of the procedure can be found in Koda 
et ah ([201 if, but here we summarize the basic concepts. The 
method consists in converting the NRO 45m map from im¬ 
age space to the uv plane, i. e., it generates the visibilities 
of the single dish map. Then, the NRO 45m visibilities are 
combined with the CARMA visibilities, and the merged uv 
sam ples are i mag ed together. Following the suggestion from 
Koda et al. ( 2011| >, we have excluded single-dish visibilities 


beyond a given mreshold in the baseline length. This thresh¬ 
old in chosen based on the the distribution of the sensitivity of 
the NRO 45m along the uv plane, and the flux recovery from 
the combined NRO 45m + CARMA map relative to the single 
dish flux calculated over regions of significant emission. For 
NGC 6946 we have cut at 7 kX , for NGC 628 we cut at 10 kX, 
and for M101 we cut at 6 kX. 

The NRO 45m CO(l —> 0) map of M101 does not cover the 
small region located in the western region of the outer disk. 
In this case, we have used the IRAM 30m CO(2 —> 1) map 
to recover t he exte nded flux using the MIRIAD task immerge 
(see Section |2.3.2| for details). 


2.3.2. Combination of CO (2 —> 1) data 

In the case of CO(2 —» 1) maps, we have used the MIRIAD 
task immerge to combine our CARMA maps with the single 
dish maps from the IRAM 30m. Immerge linearly merges a 
low resolution image with a high resolution image, which is 
equivalent to include short and large spacings in the uv plane 
into the imaging. In order to prevent this procedure to be af¬ 
fected by regions of low signal to noise ratio, we have ex¬ 
cluded edge velocity channels and we have masked the region 
where the sensitivity of the CARMA map falls below half of 
the maximum. As we did in Paper I, we have allowed im¬ 
merge to solve for the factor that put the low resolution image 
on the same flux scale as the high resolution image. In our 
combination procedure, we obtained values close to 1.05, ex¬ 
pected for data sets that have correct calibrated flux scales. 
Using regions of significant emission identified by CPROPS 
(see Section[3]for details) for both the CARMA observations, 
and the CARMA + IRAM 30m combined map, we found that 
CARMA+IRAM 30m images have typically 20% more flux 
inside the region of significant CO(2 —> 1) emission. 


2.3.3. Generation of integrated intensity maps 

We use CO spectral cubes and corresponding signal masks 
to produce integrated intensity maps. To generate a signal 
mask, we first degraded the data cube in angular and veloc¬ 
ity resolutions by a factor of 3 and 4, respectively. Then we 
identify continuous regions above 3 x tr sm0 having at least two 
pixels, where cr smo is the rms intensity of the smoothed cube. 
We further expanded each region to include any adjacent pix¬ 
els above 2 x cr smo and used the resulting regions to construct 
the signal mask. This signal mask was later used to exclude 
noise pixels in the cube, and thus generate the integrated in¬ 
tensity images. The associated uncertainty map for each im¬ 
age was derived from the error propagation through pixels in 
the signal mask. Smoothing can enhance the signal-to-noise 
ratio for extended structures in the position-position-velocity 
cube, and the dilated ma sk ap proach, similar to the one used 
in CPROPS (see Section [XT] for more details), was designed 
to help identify extended low brightness emission contiguous 
with significant emission. 


2.4. GALEXFUV data 
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RA offset (cu-csec) 


Fig. 5.— CO(2 —> 1) integrated intensity contours overlaid on the Ha: image for all the regions targeted in the northern part of NGC 628. Contours begin at 
2.5cr and are spaced by 2a. The noise at the center of the map is ~ 2 K km s _I . 
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RA offset (orcsec) 


Fig. 6.— CO(2 —> 1) integrated intensity contours overlaid on the Ha image for all the regions targeted in the southern arm of M101. Contours begin at 2.5a 
and are spaced by 2a. The noise at the center of the map is ~ 2 K km s -1 . 
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FUV data for NGC 628 and M101 are available from the 


GALEX Nearby Galaxies Survey (NGS; Gil de Paz et al. 
2007jl. The FUV band covers the wavelength range 1350 


1750 A, and the maps have an angular resolution of 5 "6. We 
applied the same procedure followed in Paper I to generate a 
sky background subtracted FUV image. We used the extinc¬ 
tion map of Schlegel et al. (1998) to correct the FUV maps 
for Galactic extinction. The FUV e xtinction w as estimated by 
using Apuv = 8.24 x E(B—V) from Wyder et al. (2007!. 


2.5. SINGS data 
2.5.1. 24 pm map 

In order to estimate the amount of star formation obscured 
by dust, we use the 24 /im data obtained with the S pitze r 
MIPS instrument provided by SINGS (Kennicutt et al. 12003) 
for NGC 628, a nd the Spitzer Local Volume Legacy (LVL, 
|Dale et al. |2009j > for M101. The angular resolution of MIPS 
is 577 (FWHM). We used the same approach we did in Paper 
I to estimate the background emission. We perform a median 
spatial filtering to remove bright spots and other features over 
a region away from the galaxy, and then measure the mean 
flux over the filtered region. 


2.5.2. Ha map 

Given the limited resolution offered by 24 pm (~ 6"). 
we need an alternative star formation tracer with a spatial 
resolution similar to that provided by our CO(2 —> 1) maps 
(~ 2"). In the case of NGC 628, an Ha image is available 
in the SINGS fourth data release as part of the SINGS an¬ 
cillary data program. The angular resolution of this Ha map 
is ~ 2", a good match to the resolution of our CO(2 —► 1) 
maps. We used the stellar continuum-subtracted Ha im¬ 
age, and we convert to flux units using the SINGS data re¬ 
lease documentation. In order to remove the [N II] contribu¬ 
tion, we have assumed ratios of [N II]A6584/Ha=0.5 and [N 
II] A6548/[N II] A6584=0.335 ( |Calzetti et al.T2005) . In the case 
of M101, we have used a Ha map observed with the Burrell- 
Schmidt telescope at Kitt Peak National Observatory (KPNO) 
( |Hoopes et al.|200T) . The map was calibrated by comparing 
the luminosities of specific H II regions to published values 
( ]Greenawalt|1998) . 

3. CLOUD PROPERTIES 


3.1. Identification 

We follow the same procedure used in Paper I to identify 
cloud structures, and we refer the reader to that paper for de¬ 
tails. The methodology is based on the cloud properties algo¬ 
rithm (CPROPS) described in |Rosolowsky & Leroy| (j2006Jl. 
CPROPS identifies clouds by isolating regions of significant 
emission that are both spatially and kinematically connected. 
More specifically, regions of significant emission are defined 
by the pixels with values greater than a threshold of n t h x o’nns 
across two consecutive velocity channels. These regions of 
significant emission are further extended to adjacent regions 
with emission greater than n e< i ge x <r rms in two adjacent veloc¬ 
ity channels. In this study we have used = 4 and w e d ge = 2 
for both CO(l —> 0) and CO(2 —► 1) maps. We keep the same 
nomenclature we used in Paper I for structures identified 
in our maps, referring to structures identified in CO(l —t 0) 
maps as “molecular complexes”, while the structures identi¬ 
fied in CO(2 —> 1) will be called “clouds”. 

The rms size o> of each identified structure is determined 
by the geometric mean of the second moments of emission 


along the major and minor axes, and the velocity dispersion 
ay is provided by the second moment of the emission along 
the velocity axis. The total flux of the structure is calculated 
by adding the flux in each pixel associated to the structure. 
A pixel with coordinates x,, y, , velocity v, and temperature 7] 
is associated to a specific structure if its position is bordered 
by the temperature isosurface defined by L e d g e, where Ledge is 
the lowest brightness temperature uniquely associated to the 
structure. We reduce the sensitivity bias by measuring these 
properties at many contour levels inside the boundary of the 
identified region and linearly extrapolating the size and the ve¬ 
locity width to the value where the brightness temperature of 
the cloud is equal to 0 K. On the other hand, the extrapolation 
applied to the flux is quadratic. The effect of finite spatial and 
spectral resolution is corrected by subtracting in quadrature 
the beam size and the spectral channel profile, approximated 
by a gaussian, from the extrapolated measurements of size 
and velocity width respectively (see Paper I for details). The 
structures identified by CPROPS with properties corrected by 
sensitivity and resolution bias are referred as resolved struc¬ 
tures henceforth in this paper. 

Following the definition of Solomon et al.' (1987]) (hence¬ 
forth S87), the cloud size is defined as R = 1.91 ay. The CO 
luminosity Leo is given by 


Leo 


Fco 


K km s 1 pc 2 K km s 1 arcsec 


180x3600 


a) 

where D is the distance to the galaxy in parsecs. The 
luminosity-based mass is obtained from Leo by using 


Mia 


M g 


: 4.4“ 


Leo 


O 


K km s pc 


y' 

2 A CO> 


( 2 ) 


where X^ Q is the assumed CO-to-fL conversion factor in units 
of 2 x 10 20 cm _2 (K kms' 1 )' 1 , and includes a factor of 1.36 


(Bell et al. 2006 Leroy et al. 2011 

Sandstrom et al. 

2013 

and theoretical (Feldmann et al. 2012J 

Narayanan et al. 

2012 


sion factor on local ISM properties such as metallicity or radi¬ 
ation field and spatial scale. Different approaches to estimat¬ 
ing the CO-to-H 2 conversion factor in normal galaxies consis¬ 
tently have found values within a factor of 2 of that found in 
the MW (~ 2 x 10 20 cm _2 (K km s -1 ) -1 ). Thus, in this study 
we assume a fixed value of Xq 0 = 1 in our mass calculations. 

Following the approach in Paper I, we calculate the virial 


mass for spheroidal clouds (see Appendix A of Bertoldi & 
McKee 1 1992 1 . Ignoring the magnetic energy and external 


pressure terms in the virial equilibrium equation, the virial 
mass is given by 



where R is the radius perpendicular to the axis of symmetry of 
the cloud. Equation |3]l corresponds to the expression for the 
virial mass assuming a spherically symmetric geometry ex¬ 
cept for a shape-dependent correction factor c/ 2 - The an factor 
is given by arcsine/e if the shape of the cloud is oblate, where 
e =(yr^) is the eccentricity of the cloud, and y = Z/R is the 
ratio between the size along and perpendicular to the axis of 
symmetry. On the other hand, If the cloud has a prolate shape. 
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then «2 = arcsinh e/e. As we did in Paper I, in the present 
work we assume that the strucutures are fully described by 
the sky-projected major and minor axis, the axis of symmetry 
is given by the major axis, and that the structures are prolate. 
We found that the mean of a 2 for CO( 1 —>• 0) complexes was 
0.8 with a rms of 0.1 for M101, and a mean of 0.8 and rms 
of 0.3 for NGC 628. On the other hand, the mean of a 2 for 
CO(2 —> 1) clouds was 0.8, with a rms of 0.1. 

We use the bootstrapping method implemented in the 
CPROPS package to estimate the uncertainties in the moment 
measurements, and this is the only source of uncertainties that 
we include in our analysis. 

3.2. Mass surface density 

We estimate the mass surface density (E H2 ) by dividing 
the luminosity-based mass by the corresponding area covered 
by the identified structure, i. e., E H2 = M\ am /(nR 2 ). Follow¬ 
ing Paper I, in order to estimate E H2 we have used the non- 
extrapolated and non-deconvolved values of the luminosities 
and areas determined at level 7) c i gc , as we intend to compare 
molecular gas surface density with the SFR tracers (24//m, 
Ha, FUV) within the boundaries of the structures identified 
by the CPROPS mask. 

3.3. Star Formation surface density 
3.3.1. SFR in CO( 1 —> 0) complexes 

Following Paper I, we estimate the star formation rate sur¬ 
face density in CO(l -+ 0) complexes using GALEX FUV and 
Spitzer 24 pm maps (see L08 for details). FUV traces the 
photospheric emission of O and B stars, measuring the star 
formation not obscured by dust over time scales of 10-100 
Myr (Calzetti et al .||2005[ l. In turn, obscured star formation is 
traced by the 24 pm emission from small dust grains heated by 
UV photons radiated from young stars. The 24 pm emission 
traces star formation over time scales of ~10 Myr (|Calzetti 
|et al.|2005) |Calzetti et al.|2007] ). L08 found that the total star 
formation surface density Esfr can be estimated by: 

-^- 2 = (8.1 X 10 _2 /f UV + 3.2 x 10- 3 / 2 4 Mm )cos/ (4) 

M 0 yr 1 kpc 

where / F uv and hA^m are in units of MJy/sr, and i is the in¬ 
clination. Equation [4] assumes continuous star formation and 
a fully sampled Initial Mass Function (IMF) over the regions 
being probed. We discuss the caveats of the se as sumptions at 
the spatial resolution of our maps in Section [53] 

In general, the amount of SF traced by 24pm dominates 
over that traced by the FUV component in the the boundaries 
of the complexes in the three galaxies. The fraction of SF 
from 24pm is 0.85 ± 0.08, 0.82 ± 0.06 and 0.81 ± 0.08 for 
NGC 6946, NGC 628, and M101 respectively. 

3.3.2. SFR in CO( 2 —> 1) clouds 

We use Ha images to estimate the amount of star formation 
in CO(2 —> 1) clouds. The Ha images offer a finer resolution 
that allows us to directly compare the star-forming regions to 
the structures identified in CO(2 —► 1) maps. As Ha emission 
can be strongly affected by dust extinction, we have corrected 
the Ha maps to account for the obscured star formation. Ac¬ 
cording to Calzetti et al. ( j2007j ), the total star formation can 
be estimated by 

-f^ r =(5.3xl0- 42 L HQ )xl0 A -/ 2 - 5 (5) 

M e yr~' 



Fig. 7.— Ssfr estimates using Ha: for CO(l —> 0) complexes as a func¬ 
tion of Ssfr predicted using a combination of FUV + 24/xm, following L08. 
Black dots represent the complexes identified in NGC 6946, while red squares 
correspond to CO(l —>■ 0) complexes detected in NGC 628 and blue triangles 
represent complexes identified in M101. Solid black lines show slopes of 2, 
1 and 0.5. We use Ahc* = 1.0 and 1.2 magnitude extinction for NGC 6946 
and NGC 626 respectively. For M101, we have applied a calibration factor to 
convert from Ha flux map to SFR using FUV + 24^/m as comparison. 

where Lu a is in units of ergs s -1 , and A Hr , is the Ha extinction. 

Due to the lack of a calibrated Ha image for M101 in the 
region we targeted, for this galaxy we have used the following 
relation to estimate the SFR 

SFR tot 

, , i — * MI a X JpVV+lAiim) \Oj 

Mq yr -1 

where F Hn is the total Ha flux in the region considered, and 
/FUV+ 24 /jm is a calibration factor to convert from the Ha map 
units to SFR rate using the FUV+24/itn map as a reference 
for the total star formation rate. 

If we assume that the SFR values estimated using FUV+ 
24pm are correct, we can estimate the Ha extinction Ah q for 
NGC 628 and the calibration factor /fuv+ 24 /mii for M101 by 
comparing the star formation rate calculated using Equations 
((5) and Q to the values derived from Equation Q. In both 
cases, the comparison was performed by eye. Figure [7] shows 
the comparison between the different methods to estimate the 
SFR. For NGC 628, we found that a mean value of A Ho = 1 .2 
applied to Equation © roughly recovers the SFR traced by 
FUV+ 24pm. 

By assuming a constant value of An a for NGC 6946 and 
NGC 628, and a fixed calibration factor for M101, we may 
be introducing some scatter in our SFR measurements for 
CO(2 —► 1) clouds. According to Figure IT] there is a scatter 
of a factor of ~ 2 about the identity line between Esfr from 
Ha and the Esfr from FUV+24// m. This scatter is similar to 
that shown in Figure 23 of L08. However, we notice that the 
Esfr from FUV+24 p m consistently overestimates the Esfr 
from Ha for Esfr 0.01 M 0 yr _1 kpc -2 . We further discuss 
this discrepancy in the context of a cirrus component in 24pm 
in Appendix [A] 

4. RESULTS 

4.1. Statistics of identified structures 

The number of identified structures in each galaxy by 
CPROPS is shown in Table [3] along with the number of 
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TABLE 3 

Mean values of resolved structure properties. 


Property 

NGC 6946 

NGC 628 

M101 

CO( 1 — y 0) complexes 

No. identified structures 

43 

34 

35 

No. resolved structures 3 

33 

19 

17 

{R)( pc) 

180 ±9 

142 ± 12 

162 ± 12 

(<t v ) (km s -1 ) 

6.7 ± 0.3 

4.8 ± 0.4 

5.7 ±0.4 

(Leo) (10 5 K km s _1 pc 2 ) 

20 ± 1 

11 ± 1 

15 ± 1 

(M vir >(lO 6 M 0 ) 

11 ±2 

5 ± 1 

7 ± 1 

(S H 2> (Mq pc -2 ) 

44 ±2 

36 ±2 

49 ±3 

(S S fr)( 1O‘ 3 M 0 yr 1 kpc' 2 ) 

36 ±2 

21 ±2 

62 ±7 

CO(2 — >■ 1) clouds 

No. identified structures 

64 

28 

52 

No. resolved structures 3 

34 

25 

38 

(if) (pc) 

65 ±4 

98 ±6 

96 ±6 

(cr v ) (km s -1 ) 

5.7 ±0.3 

5.1 ± 0.3 

5.4 ±0.3 

(Leo) (10 4 K km s -1 pc 2 ) 

40 ±2 

59 ±3 

46 ±3 

(M vir )(lO 5 M 0 ) 

24 ± 1 

33 ±4 

37 ±5 

(S H 2> (Mq pc' 2 ) 

70 ±3 

52 ±3 

37 ± 2 

(Esfr)(1O' 3 M 0 yr' 1 kpc' 2 ) 

90 ±6 

62 ±5 

134 ± 15 


Notes. 


a Structures with properties corrected by sensitivity and resolution bias. 


structures with measured values of R and <r v after correc¬ 
tion by sensitivity and resolution bias is applied (i. e., re¬ 
solved structures). Also, Table [3] shows the mean values of 
the properties for resolved structures. The CO(l —> 0) com¬ 
plexes have resolved radii between 23 pc and 426 pc, and ve¬ 
locity dispersions between 2.6 km s' 1 and 13.1 km s' 4 . The 
luminosity Leo for complexes has values between 1.3 xlO 5 
K km s' 1 pc 2 and 108 x 10 5 K km s 1 pc 2 , and virial masses 
from 4.7 x 10 5 M 0 to 766 x 10 s Mq. 

On the other hand CO(2 —y 1) clouds have resolved radii 
between 25 pc to 171 pc, and velocity dispersions between 
1.2 km s _I to 9.6 km s -1 . Leo for clouds has values between 
2.7 x 10 4 K km s -1 pc 2 and 280 x 10 4 K km s -1 pc 2 , and virial 
masses from 7.6 x 10 4 M 0 to 1810 x 10 4 Mq. 

At scales sampled by the CO(l —>• 0) observations (~ 150- 
200 pc), the mean surface density for each galaxy, (Shi), is 
in the range of ~ 35-50 M 0 pc' 2 . In contrast, the higher 
resolution CO(2 -A 1) observations reveal some differences in 
the mean values. NGC 6946 shows the highest mean of ~ 
70 M & pc' 2 for the CO(2 — y 1) clouds followed by NGC 628 
with a mean of ~ 50 M 0 pc“ 2 , and finally M101 shows the 
smallest mean surface density ~ 40 M 0 pc 2 . This difference 
in the (Sea) of NGC 6946 basically reflects the fact that the 
overall CO brightness of NGC 6946 is higher than the other 
two galaxies, M101 and NGC 628. 

Overall, the Esfr values for the regions observed in the 
three galaxies cover similar ranges, from log( ^ Ss f k - 2 ) ~ 

-2.7 to 0 for both complexes and clouds. We identify 
some complexes in NGC 6946 and M101 with Esfr > 
0.1 Mq yr' 1 kpc 2 . These complexes are located in on-arm re¬ 
gions, where local star formation activity is higher compared 
to other areas surveyed over the three galaxies. 

At higher resolution, clouds follow a behavior simi¬ 
lar to the complexes, with clouds located in specific re¬ 
gions of the spiral arms exhibiting active star formation (> 
0.1 Mq yr kpc ") constituting the high end of the Esfr val¬ 
ues. 


4.2. Scaling relations 
4.2.1. Bayesian inference method 

We have used a Bayesian inference method to find the best 
fit relations between the properties of the structures. This ap¬ 
proach differs from the bisector linear regression method im¬ 
plemented in Paper I for NGC 6946. Bayesian fitting meth¬ 
ods have been successfully employed in several astrophysical 
analyses, such as the derivation of the extinction law in the 
Perseus molecular cloud complex Foster et al.|2013| l. Type la 
supernova light curve inference (Mandel et al.|2011]), the size¬ 
line width relati on in the de nse ISM of the Central Molecu¬ 
lar Zone ( Shetty et al.|2012| , and the extragalactic Kennicutt- 
Schmidt rela tion ( |Shetty et al. 20131. We refer the reader to 
|Kelly| ( [2007| l for details of the Bayesian regression fit method, 
and we limit our description here to the basic concepts. The 
Bayesian approach estimates the joint posterior probability 
distribution of the regression parameters given the observed 
data, and draws the error in each measured quantity from 
some a priori defined distribution which should reflect the un¬ 
certainties in the measurements. Bayes’ theorem yields the 
posterior distribution of the parameters 9 given the observed 
data (x,y) as 


p(0\x,y) oc p(x,y\6)p(6), 


(7) 


where p{0) is the prior distribution of the parameters, and 
p(x,y\9) is the probability of the data given the parameters 
6. Although a direct derivation of the posterior probability 
distribution is computationally expensive, the Markov chain 
Monte Carlo (MCMC) routine allows one to sample the prob¬ 
ability distribution of the fit parameters through random draws 
in parameter space. This yields a histogram of the marginal 
probability distribution, allowing the estimation of the me¬ 
dian and error for each parameter. In our analysis, we have 
used the Gibbs sampling method for generating random draws 
from the posterior distribution ( Kelly [2007} . As this method 
accounts for the uncertainties of the dependent and indepen¬ 
dent variables of the fit, it provides more realistic uncertainty 
estimates for each parameter of the regression. 

Through this paper, we model the scaling relations using the 
standard power law form given by y = ax a . Following the pre¬ 
scription introduced by |Kelly| (J2007J), the fitting is performed 
in the log space. 


log (y) = A + a log (x) + e scat , 


( 8 ) 


where A = log (a), and e sca t is the scatter about the regression 
line. e scat is assumed to have mean of 0 and dispersion of < 7 . 
Thus, the three parameters involved in the fitting method are 
A,a and o. For each scaling relation, we have run 2 x 10 4 
random draws to sample the probability distribution of the fit 
parameters. The Bayesian regression fit method used in this 
paper (jKelly||2007|) allows us to include the non-detections in 
the y-axis only. However, in Section 4.3.3 we will introduce a 
new approach to include Eh 2 non-detections in the context of 
the K-S relation. 


4.2.2. Size-Line width relation 

Figure [8] shows the scatter plot of the line width and size 
of CO(l —► 0) emitting complexes on a logarithmic scale. For 
comparison, we include the relations found by S87 and B08 in 
the figure. Because we are interested in the variations of the 
regression fit parameters at size scales representative of the 
values observed in complexes and clouds, we have normalized 
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TABLE 4 

Bayesian regression parameters for CO(l —>■ 0) complex scaling relations. 


Relation 

a 

90% HDI 

A 

90% HDI 

(7 

90% HDI 

NGC 6946 

(7 v R 

0.23 

[-0.13,0.67] 

0.79 

[0.61,0.88] 

0.09 

[0.06, 0.14] 

Mvu-Lco 

1.15 

[0.84, 1.40] 

6.56 

[6.38, 6.71] 

0.16 

[0.09, 0.25] 

NGC 628 

<r v —R 

0.05 

[-0.35, 0.50] 

0.70 

[0.57, 0.78] 

0.06 

[0.01,0.12] 

Mvh-Lco 

1.17 

[0.64, 1.63] 

6.42 

[6.21, 6.54] 

0.09 

[0.03, 0.25] 

M101 

(t v —R 

0.19 

[-0.20, 0.61] 

0.70 

[0.57, 0.81] 

0.08 

[0.03,0.17] 

Mvu-Lco 

1.53 

[1.16, 2.18] 

6.17 

[5.86, 6.35] 

0.11 

[0.04, 0.32] 





Fig. 8.— Size-line width relation for the CO(l —»■ 0) complexes found in this study (top-left). Black dots represent the complexes identified in NGC 6946, 
while red squares correspond to complexes detected in NGC 628 and blue triangles represent complexes identified in M101. The dashed line illustrates the 
relation found by S87, and the dashed dotted line represents the fit found by B08. The marginal distribution of the slope a from the Bayesian regression analysis 
is shown in the top-right panel, the intercept coefficient A is shown in the bottom-left panel and the dispersion a is shown in the bottom-right panel, with 

lo S( A + “ lo §( TOiTpc ) + ■ 


the size variable by an appropriate value (100 pc). Thus, for 
the size-line width relation, we have expressed the sizes in 
units of 100 pc, so the linear regression is performed over the 
relation 

108 (k^Fr) =A+al08 (io^) +e “- <9) 

The distributions of the regression fit parameters using the 
Bayesian inference approach are shown in Figure [HI The dis¬ 
tributions of the slope a, the y-axis intercept coefficient A, 
and the dispersion of the scatter a for the size-line width rela¬ 
tions show overall consistent behaviours for the three galaxies 
in our sample. The peak of the slope distribution is < 0.5 for 
the three galaxies, and the wide range of the HDI (High Den¬ 


sity Interval, defined as the interval that encloses the 90% of 
the probability distribution) for the slopes (HDI ~ 0.8) indi¬ 
cates the a parameter is poorly constrained, i. e., we do not 
find a strong correlation between size and line width for the 
complexes. For NGC 628, the complexes show line widths 
smaller than 7 km s _l for sizes that span a range of 50 to 300 
pc. Among the three galaxies in our sample, NGC 628 shows 
the smallest inclination (7°), which suggests that the CO ve¬ 
locity dispersion observed in this galaxy is less affected by 
cloud blending or superposition of several clouds along the 
line of sight. The intrinsic scatter about the linear relation 
is similar for the three galaxies, although the distribution of 
er is shifted to smaller values for NGC 628. The values of 
the peaks and the HDI provided by the Bayesian inference 
method for the slope a, the intercept coefficient A, and the 
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TABLE 5 

Bayesian regression parameters for CO(2 —> 1) clouds scaling relations. 


Relation 

a 

90% HDI 

A 

90% HDI 

<J 

90% HDI 

NGC 6946 

u v —R 

0.27 

[-0.11,0.67] 

0.83 

[0.75, 0.87] 

0.09 

[0.06, 0.14] 

o 

i 

1.17 

[0.94, 1.46] 

6.64 

[6.53, 6.74] 

0.19 

[0.11,0.31] 

NGC 628 

<t v —R 

-0.17 

[-0.77, 0.47] 

0.73 

[0.67, 0.76] 

0.06 

[0.01,0.10] 

M v ii~Lco 

0.65 

[0.34, 1.06] 

6.66 

[6.55, 6.74] 

0.12 

[0.04, 0.22] 

M101 

<t v —R 

0.37 

[-1.48, 2.78] 

0.76 

[0.67, 0.80] 

0.08 
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Fig. 9.— Size-line width relation for the CO(2 -A 1) clouds found in this study (top-left). Symbols and lines are the same as FigureJS] The panels showing the 
marginal distribution of the slope a, the intercept coefficient A and the dispersion of the scatter of the relation a from the Bayesian regression analysis are placed 
in the same way as Figurejs] 


dispersion <j for complexes are shown in Table HI 
Figure [9] shows the a v - R plot for the CO(2-% 1) clouds 
found in our sample of galaxies, and the values for the lin¬ 
ear fit using Equation [9] are shown in Table [5] Similarly to 
the parameter distribution for complexes shown in Figure [8] 
for clouds we observe that the distribution of the three pa¬ 
rameters overlaps substantially for the three galaxies in our 
sample. Although the slope distribution peaks at similar val¬ 
ues in both NGC 6946 and M101 (~ 0.32), the HDI for p(a) 
for M101 (~ 2.8) is wider than NGC 6946 (~ 0.8). The inter¬ 
cept coefficient A is larger for NGC 6946 clouds than for the 
M101 counterparts, indicating a larger mean velocity disper¬ 
sion for clouds with similar sizes located in the NGC 6946. 
Similarly to the behaviour observed at larger size scales in 
complexes, NGC 628 shows a nearly flat distribution of the 
size-line width relation (a pe ak ~ -0.1 ) for clouds. The line 
widths are smaller than 7 km s -1 , and generally smaller than 


the velocity dispersions observed in clouds with similar sizes 
located in the other two galaxies. 


4.2.3. Virial mass-CO Luminosity relation 

The virial mass-luminosity relations for the complexes lo¬ 
cated in the three galaxies are shown in Figure[T0] We observe 
that the complexes roughly fall along the relation defined by 
a constant CO-to-FE factor Xco = 2 x 10 20 cm _2 (K km s -1 )~', 
which was adopted in Equation d2l). Nevertheless, we no¬ 
tice that the complexes in M101 and NGC 628 show slightly 
smaller virial masses than NGC 6946 for the same CO lu¬ 
minosities. Assuming virial equilibrium at such large scales, 
this difference can be interpreted as smaller values of Xco in 
the regions surveyed in NGC 628 and M101 than the bright 
eastern part of NGC 6946, This is con sistent with the val¬ 
ues reported by Sandstrom et al. (2013) using dust maps as 
the tracer of the mass. For the same regions we observed 
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Fig. 10.— The Virial mass-Luminosity relation for the complexes found in this study (top-left). Solid blue lines show different AY o values. Down arrows 
illustrates upper limits of unresolved clouds. The marginal distribution of the slope a from the Bayesian regression analysis is shown in the top-right panel, the 
intercept coefficient A is shown in the bottom-left panel and the dispersion a is shown in the bottom-right panel, with log( ) = A + alog( |(|6 ) + £scat- 


in NGC 628 and M101, they found mean values of the CO- 
to-FF conversion factor slightly below (Xco ~ 1.4 - 1.7 x 
10 20 cm -2 (K km s -1 ) -1 ) the value used here. 

Similarly to the size-line width relation, we have normal¬ 
ized the CO luminosity of the complexes by a representative 
observed value, chosen to be 10 6 K km s _1 pc 2 , to find the best 
fit for the virial mass-luminosity relation . Then, the relation 
utilized for the Bayesian regression fit is given by 


, /M vir 

losl W 


-A + a log 


Leo 


10 6 Kkm s _1 pc 2 


+ c sca t, (10) 


where Leo is the CO( 1 —> 0) luminosity for complexes, and 
the CO(2 —> 1) luminosity for clouds. 

Figure [TO] shows the distributions of the linear fit param¬ 
eters using the Bayesian inference method. Following the 
procedure detailed by |Kelly ( f2007j l, we have incorporated the 
upper limits of the virial masses in the fitting from spatially 
unresolved clouds. According to this procedure, an additional 
indicator is introduced in the statistical model, which indi¬ 
cates whether a data point in the dependent variable is cen¬ 
sored (i. e., non-detected) or not. The fit is close to a lin¬ 
ear relation for galaxies NGC 6946 and NGC 628, but with 
the latter showing more scatter in the slope distribution likely 
due to the small range in velocity dispersion observed in this 
galaxy (see Figure[8|. M101 shows a slope peak a pea k ~ 1.5, 
and the A pea k intercept coefficient is smaller than the other two 
galaxies. Moreover, NGC 628 shows a smaller A coefficient 
than NGC 6946, consistent with the smaller Aco inferred for 
NGC 628 and M101. This distribution of the intrinsic scatter 


about the relation is similar for the three galaxies. We have 
compared the values found in this paper to the values found 
by previous works. The complexes identified in NGC 6946 
are located along the —Leo relation found by B08, but 
are systematically above the relation found by S87 for Galac¬ 
tic clouds. On the other hand, the low luminosity complexes 
identified in NGC 628 and M101 are closer to the Galactic 
values found by S87, while high luminosity complexes tend 
to be closer to the extragalactic values found by B08. 

The virial mass-luminosity relation for the CO(2 —► 1) 
clouds is shown in Figure [IT] The clouds are roughly consis¬ 
tent with the CO-to-FF conversion factor used here, and con¬ 
sistent with the relation found by B08 for extragalactic clouds. 
As for complexes, here we have incorporated the upper limits 
in the fitting from unresolved clouds. The slope distributions 
are slightly different for the three galaxies, with M101 show¬ 
ing the steepest relation (a pea k ~ 1.5) and NGC 628 showing 
the most shallow one (a pea k ~ 0.6). 


4.3. The Star Formation Law 

In this section we show the relation between the molecu¬ 
lar gas and star formation surface density in the regions we 
observed in NGC 6946, NGC 628 and M101. We begin by 
showing the resulting relations using the boundaries of the 
CO(l —>• 0) complexes and CO(2 —> 1) clouds to define the 
area over which we calculate the quantities. Then, we com¬ 
pare these results with the K-S relation calculated by using a 
uniform grid across the regions observed in the galaxies in our 
sample. 
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Fig. 11.— The Virial mass-Luminosity relation for clouds (top-left). Solid blue lines show different Xco values. Down arrows illustrates upper limits from 
unresolved clouds. The panels showing the marginal distribution of th e slo pe a, the intercept coefficient A and the dispersion of the scatter of the relation a from 
the Bayesian regression analysis are placed in the same way as FigureflO] 


4.3.1. Eh2 vs. Esfr relation based on identified structures 

The left panel of Figure [l2| shows the S H 2 vs. £$fr relation 
for the CO( 1 —>• 0) complexes identified in the three galax¬ 
ies. We observe that the complexes are, in general, located in 
the same area of the £ H 2 vs. £$fr diagram. However, com¬ 
plexes located in NGC 6946 seem to show a stronger Eh 2 vs. 
Esfr correlation than the structures identified in the other two 
galaxies. Additionally, there is a subset of complexes in M101 
with SFR several times higher than the SFR values found for 
similar molecular gas surface densities found in NGC 6946 
and NGC 628. 

The Bayesian regression fitting to the £h2 vs. Esfr relation 
is performed similarly to the fit performed on the scaling re¬ 
lations presented above. In this case, we have normalized the 
molecular gas surface density by the representative value of 
50 Mg pc~ 2 . Thus, the relation used for the fitting is given by 


log 


Esfr 


M q yr“' kpc 


-2 


= A + a log 


£h2 


50 Mq pc 


-2 


3“^scat- 


( 11 ) 

The distributions of the slope a , the intercept coefficient A 
and the dispersion of the intrinsic scatter a of the regression fit 
for CO(l -4 0) complexes are displayed in Figure[l2] Table[6] 
shows the values of the regression parameters for complexes. 
For NGC 6946, the distribution of the slope has a maximum 
at 1.33. On the other hand, the peaks of the slope distribu¬ 
tions are ~ 1.4 and 2.5 for NGC 628 and M101 respectively. 
In the case of the NGC 6946, the HDI of the distribution of 
a covers a range of ~ 0.7, extending from 0.97 to 1.65. For 


NGC 628, the HDI covers a range of ~ 3.5, extending from 
0.43 to 3.95. Finally, the distribution of the slope of the K-S 
relation for complexes in M101 shows the largest range ex¬ 
tending from 1.09 to 6.77. Thus, the distribution of the slope 
is consistent with a super-linear relation for M101. In the case 
of NGC 6946 and NGC 628, although the largest area of the 
distribution of the slope is located above 1, the HDI extends 
to sub-linear values, especially for NGC 628. By comparing 
the distributions of the intrinsic scatter about the Eh 2 vs. Esfr 
relation for the three galaxies, we observe that NGC 6946 of¬ 
fers the tightest relation among the three galaxies in the sam¬ 
ple (peak of er ~ 0.17), although significant overlap with the 
distribution of the intrinsic scatter in NGC 628 is seen. On 
the other hand, the intrinsic scatter is larger in M101 (peak 
of a ~ 0.45) due to the large range of SFR values observed 
in CO(l -4 0) complexes, in contrast to the narrow dynamic 
range in molecular gas surface density values. 

Figurefl3lshows the relation of E H 2 vs. Esfr for CO(2 -4-1) 
clouds. In Table [6] we show the values of the regression 
parameters for clouds. The general trend is similar to that 
observed for complexes: a tighter correlation holds for the 
clouds identified in NGC 6946. The distributions of the slope 
show peaks at 1.57, 1.98 and 2.58 for NGC 6946, NGC 628 
and M101 respectively. Overall, clouds show similar E H 2 vs. 
Esfr relations to those found for large structures identified as 
complexes. In the case of NGC 6946 and M101, the HDI of 
the probability distribution of slopes is consistent with super- 
linear relations. For NGC 628, the HDI extends to values be¬ 
low 1. By comparing the intrinsic scatter for complexes and 
clouds shown in Figure [12] and [13] respectively, we observe 
that the scatter about the Eh 2 vs. Esfr relation increases for 
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TABLE 6 

Bayesian regression parameters for both CO(l —* 0) complexes and CO(2 —» 1) clouds E H 2-£sfr relation. 


Relation 

a 

90% HDI 

A 

90% HDI 

a 

90% HDI 

NGC 6946 

Eh 2 — Esfr CO(l —> 0) complexes 

1.33 

[0.97, 1.65] 

-1.48 

[-1.54,-1.41] 

0.17 

[0.11,0.23] 

Spi 2 — Esfr CO(2 —y 1) clouds 

1.57 

[1.06,2.02] 

-1.51 

[-1.63,-1.39] 

0.38 

[0.31,0.47] 

NGC 628 

Sh 2 _ Esfr CO(l —» 0) complexes 

1.38 

[0.43, 3.95] 

-1.59 

[-1.78,-1.26] 

0.24 

[0.13,0.34] 

Spj 2 — ^sfr CO(2 —y 1) clouds 

1.98 

[0.96, 3.09] 

-1.54 

[-1.75,-1.37] 

0.51 

[0.39, 0.70] 

M101 

Eh 2 — ^sfr CO(l —> 0) complexes 

2.48 

[1.09,6.77] 

-1.65 

[-1.84,-1.47] 

0.45 

[0.19, 0.59] 

Epi 2 — Esfr CO(2 —y 1) clouds 

2.58 

[1.81,3.71] 

-1.13 

[-1.34. -0.94] 

0.58 

[0.44, 0.73] 






Fig. 12.— Star formation rate vs. molecular gas surface density relation for CO(l —>■ 0) complexes (top-left). The SFR is calculated from the FUV + 24/zm 
maps. Black dots represent the complexes identified in NGC 6946, red squares correspond to CO(l —»■ 0) complexes detected in NGC 628 and blue triangles 
represent complexes identified in M101. The vertical dashed line illustrates the sensitivity limit of our maps, ~10 M@ pc 2 . Black dot-dashed lines show 
constant molecular gas depletion times (SFE -1 ) of 0.1, 1 and 10 Gyr. The marginal distribution of the slope a from the Bayesian regression analysis is shown 
in the top-right panel, the intercept coefficient A is shown in the bottom-left panel, and the dispersion of the scatter a is shown in the bottom-right panel, with 


log( 


^SFR 


A/q yr 1 kpc 


i)=A + a log( 


50 Mq pc 


r) + Cscat- 


CO(2 —>• 1) clouds. More scatter could have been introduced 
by assuming a constant Ha extinction factor in our SFR cal¬ 
culations for CO(2 —>• 1) clouds as was mentioned in Section 
3.3.2| Additionally, the higher scatter in the K-S relation for 
CO(2 -A 1) clouds may be the result of a relatively narrow 
range of observed CO(2 —1 1) brightness, added to the intrin¬ 
sic scatter introduced to the relation by local variations of the 
molecular gas surface density and star formation at smaller 
scales than those traced by CO(l —> 0) observations. Given 
the small regions sampled by our COG —> 0) and CO(2 —1 1) 
observations (with a bias to regions of bright CO emission), it 
is not possible to generalize the Eh 2 vs. E$fr relation found 
for complexes and clouds as a global trend. In the next section 
we will compare the relations that we have found based on 


CO( 1 —> 0) complex boundaries with the relation derived by 
using a uniform grid covering the COG —>■ 0) and SFR maps, 
and we investigate the effect of non-detections in the x-axis of 
the K-S relation on the derived linear relations. 

4.3.2. £h 2 vs. Esfr using a uniform grid 

So far, we have used the sizes of the identified structures 
(complexes and clouds) to define the area over which to cal¬ 
culate the properties. That approach is desired (and needed) 
when the properties require the knowledge of the size and line 
width of the region under study. Nevertheless, comparison of 
star formation and molecular gas surface density require some 
other considerations. For instance, the Eh2 vs. Esfr relation 
is affected by sensitivity bias in both axes. In our CO(l —> 0) 
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Fig. 13.— Star formation rate vs. molecular gas surface density relation for CO(2 —> 1) clouds (top-left). Symbols and lines are the same as in Figure |l2] 
except the sensitivity limit is ~15 Mq pc 2 . The marginal distribution of the slope a, t he intercept coefficient A, and the dispersion of the scatter cr from the linear 


Bayesian regression analysis are shown in the same panel configuration as Figure |12j with log( 


S SFR 


v Mq yr 1 kpc' 


^)=A + a log( 


50 Mq pc 


r) + €scat- 


maps, our Eh 2 detection limit is ~ 10 Mq pc ~ 2 , so in our 
survey we may be missing low-brightness molecular clouds, 
more likely to be located in the inter-arm regions. Several sta¬ 
tistical approaches have been implemented to correct for such 
bias ( |Leroy et al.|2013| |Blanc et al,|2009) . Additionally, in- 
cluding the full coverage of the region with observations is an 
important factor in the observed E H 2 vs. Esfr relation. Defin¬ 
ing sample regio ns using the stmcture boundaries (as we did 
in Section [4.3.1 1 may lead us to miss of some CO emission 
that is not attributed to any identified complex or cloud. An 
additional bias could be introduced by selecting young clouds 
with incipient star formation and abundant molecular gas, but 
excluding more evolved star-forming regions that have de¬ 
stroyed their natal clouds and are weak in CO emission. 

The common approach to study the E H 2 vs. Esfr relation 
in galaxies has been t he use of a uni form sampling grid across 
the region observed (Bigiel et al.||2008| |Leroy et al.||2013[ 
|Blanc et al.|2009). This approach has the advantage of cover¬ 
ing uniformly the region under study, allowing the inclusion 
of CO non-detections into the analysis. In order to compare 
the Eh2 vs. Esfr relation for complexes with the relation ob¬ 
tained by using a uniform grid across the regions observed, 
we have generated gridded maps of SFR and molecular gas 
surface density. The Esfr ma P is created using Equation 
Q. i. e„ we use the combination of FUV and 24p.m. On 
the other hand, the Eh 2 map is derived from the CO( 1 — >• 0) 
integrated intensity map. We have sampled the maps using 
a hexagonal grid such that each surface element considered 
in the analysis is a hexagon with edge size of 8", so each 
grid elements contains ~ 1.5 - 2 times the beam size of the 


maps. This edge size corresponds to 213 pc for NGC 6946, 
283 pc for NGC 628 and 287 pc for M101. Every individ¬ 
ual grid point in the sample is the average of the pixel values 
inside the corresponding hexagon. The re-gridded maps of 
Eh 2 and Esfr for NGC 6946, NGC 628 and the two regions 
observed in M101 are shown in Figure [T4| Figure |1~5| shows 
the Eh 2 vs. Esfr scatter plot for the molecular complexes and 
for the uniformly sampled grid. We observe that the com¬ 
plexes follow the distribution of the grid values for the three 
galaxies at high values of Eh 2 - The uniform grid has intro¬ 
duced a group of data points with Eh 2 in the range ~ 10-50 
Mq pc' 2 , but with a narrow range in Esf r lik ely to be asso¬ 
ciated with the inter-arm zone. In Section [5741 we discuss the 
differences in the star formation properties of the clouds in 
the context of on-arm and inter-arm regions. Figure p~5] also 
shows the distribution in the SFR axis of Eh 2 non-detections 
for each galaxy. Assuming a line width window of ~ 7 km s' 1 
to recover the entire flux, we have estimated our Eh 2 non¬ 
detection threshold to be 10 M 0 pc' 2 . We have used the same 
value for all the galaxies in our sample given that the sensitiv¬ 
ity is similar across our maps (~ 0.4 K). We notice that while 
in NGC 628 the Eh 2 non-detections show SFR values below 
~ 0.03 Mq yr' 1 kpc' 2 , Eh 2 non-detections can reach higher 
SFR values of up to ~ 0.2 Mq yr' 1 kpc' 2 for the galaxies 
NGC 6946 and M101. In the next section, we investigate the 
effect of non-detections on our Bayesian regression analysis 
of the K-S relationship. 


4.3.3. Eh 2 non-detections 
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Fig. 14.— Left: Molecular gas surface density for CO(l —> 0) map using a uniform grid across the the areas observed in NGC 6946 (top), NGC 628 (middle) 
and M101 (bottom). Color bar is in units of Mq pc~. Right: Star formation surface density using the same uniform grid approach as in the right panel for NGC 
6946 (top), NGC 628 (middle) and M101 (bottom). Color bar is in units of Mq yr 1 kpc 2 . 


The common approach to dealing with non-detections in 
previous studies of the K-S relationship has been to clip the 
data below a given sensitivity threshold in the CO flux. How¬ 
ever, as such clipping can lead to a significant bias in fitting 
a line to the log(EH 2 )-log(EsFR) relationship, recent studies 
have proposed different ways to inc lude the non-detections in 
the K-S analysis (for example, see ILeroy et al.||2013| |Blanc 
let al.|2009). 


In this paper, we have included the molecular gas surface 
density non-detections in the Bayesian analysis by assum¬ 
ing that Eh 2 grid values for each galaxy follow a log-normal 
distribution at a given Esfr- We adopt this model largely 
for its simplicity, but note that log-normal density distribu¬ 
tions arise naturally in simulations of interstellar turbulence 
(e. g.,|Wada & Norman|2001jl, although the width of the log- 
normal distribution may vary with the Mach number (e. g., 


| Passot & Vazquez-Semadeni|1998| l and/or star formation rate 
( |Kravtsov||2003| ). Because we do not know the Eh 2 values 
for non-detections, we can only use the distribution of de¬ 
tected E H 2 values to infer the underlying log-normal distribu¬ 
tion. We divide the distribution of observed Esfr into a series 
of bins in log(EsFR)- In order to assure a significant num¬ 
ber of Eh 2 detections in each bin, we have selected a width of 
Alog(EsFR) = 0.5. We allow a 50% overlapping between bins 
to have a smooth transition of the E H 2 log-normal distribution 
parameters across the range of Esfr values. Our approach is 
only applied to log(EsFR) bins for which the fraction of Eh 2 
detections is larger than 10%. 

For each log(EsFR) bin, we fit the erfc function to the frac¬ 
tional number of grids with detected molecular gas surface 
density larger than a given log(E H 2 ) value. 
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Fig. 15. — Eh 2 vs. Esfr relations using a uniform grid (black dots), and using the CO(l —> 0) complexes (red squares) for the three galaxies in our sample 
NGC 6946 (left), NGC 628 (middle) and M101 (right). Black arrows illustrate the £h 2 non-detections at the sensitivity limit of 10 Mq pc -2 . 


/(> log(E H2 )) = 


N(> log(£ H 2)) 

Mo, 


( 12 ) 


where N(> log(E H 2 )j is the cumulative number of grid points 
with molecular surface density larger than log(E H 2 ), and Mot 
is the total number of grid values in the sample. The values 
of the fractional distribution /(> log(E H2 )) and their uncer¬ 
tainties are estimated in two steps. The first task is to estimate 
the adopted value and its uncertainty for each bin in the dif¬ 
ferential M(log(£ H 2 )) distribution. By assuming that the un¬ 
certainties in the Eh 2 detections are Gaussian distributed, we 
generate a new set of observed values Eh 2 using 


- £h2 + ^£ 


J H2 


03) 

where Eh 2 is the observed value, and ,)E H 2 is drawn from 
a normal distribution with mean 0 and dispersion equal to 
the error of the observed Eh 2 . For each generated sample 
of Eh 2 values, a new histogram M(log(En 2 )) is built with a 
fixed bin size of 0.1. The adopted M(log(En 2 )) value for each 
bin and its uncertainty is calculated from the mean and stan¬ 
dard deviation respectively, over 1000 different realisations of 
M(log(E H2 )). 

The second step is to estimate the values and uncertain¬ 
ties of each bin in the fractional cumulative distribution /(> 
log(En 2 ))- They are calculated by generating 1000 different 
realisations of /(> log(En 2 ))- In this case, each realisation 
of /(> log(En 2 )) is the result of applying Equation [13] to the 
adopted value of M(log(E H2 )) and derived uncertainty. 

Finally, we find the best fit erfc function to the /(> 
log(E H2 )) distribution for each log(EsFR) bin. To perform the 
fit, we have used the MPFIT fitting package given by Mark- 
|wardt| ( j2009j ). The non-detections are randomly distributed 
using a log-normal distribution in Eh 2 with mean and disper¬ 
sion consistent with the erfc fit for each Esfr bin. As an exam¬ 
ple of our methodology, Figure[l 6 ]shows the cumulative func¬ 
tions and complete distribution of both Eh 2 non-detections 
and detections for each log(EsFR) bin in the galaxy NGC 628. 
In Figure [TT] we show the Eh 2 vs. Esfr scatter plot for the 
three galaxies including both detections and non-detections 
of E h2 . 

Once we have assigned values to Eh 2 non-detections, we 
can include them i n the B ayesian analysis of the K-S relation¬ 
ship. As in Section 4.3.1 we have used 2 x 10 4 random draws 
to derive the posterior probability distribution of the linear 
fit parameters in the Bayesian approach. This procedure was 


performed for the following cases: considering only E H 2 de¬ 
tections and including detections + non-detections. In Figure 
T9|we show the probability distributions of the parameters a, 
4and a from the Bayesian regression fit for each case and for 
each galaxy in our sample. Table[7]shows the peak and 90 % 
HDI of the probability distributions for each fitting parameter. 
Representative relations built by using the peak values of each 
parameter distribution are shown in Figure [17] 

In the case of NGC 6946, the slope a distribution peaks at 
0.63 when only Eh 2 detections are considered in the fit. On 
the other hand, if all the non-detections are included in the fit, 
the slope distribution peaks at 0.38. The distributions of the 
slope for both cases do not show a significant overlapping. 
On the other hand, the distributions of the coefficient A and 
the dispersion of the intrinsic scatter a when only detections 
are considered and when both detections and non-detections 
are considered are consistent. 

In contrast to NGC 6946, the distribution of the slope a 
for NGC 628 does not show a significant change when the 
Eh 2 non-detections are included or not to the fitting of the K- 
S relation. In both cases, the peak of the distribution of the 
slope is 0.63, and the 90 % HDI is slightly larger in the case 
when only detections are considered. The distribution of the 
coefficient A is consistent between the two cases. 

For M101, the distribution of the slope peaks at 0.98 when 
only Eh 2 detections are considered in the fit, while when the 
non-detections are included the peak is 0.83. However, both 
distributions overlap significantly, mainly due to the broad 
distribution of slopes in the detection-only case. Further evi¬ 
dence of the poor correlation in the K-S relation for M101 is 
the large scatter a yielded by our fitting approach when only 
detections are considered in the analysis. 

4.3.4. Cirrus component of the Esfr 

So far, our analysis has not considered the systematic un¬ 
certainties in the Esfr measurements. In their analysis of of 
the IRAM HERACFES CO survey, [Leroy et ai.’j (j2012j) dis¬ 
cuss the bias inherent to each component of the commonly 
used “hybrid” tracers FUV+24//m and Ha+ 24/<m. Addition¬ 
ally, they estimate the contamination of the 24 pm band by 
emission not associated with current star formation by using 
physically motivated dust models. At 1 kpc resolution, they 
find that systemic differences among these different tracers 
are less than a factor of two for Esfr > 1 O - 3 M 0 yr 1 kpc -2 . 

Recent works on the extragalactic K-S law have determined 
that the index in the power relation between molecular gas 
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Fig. 16.— Left column: Cumulative distributions, /(> log(En 2 )) ? for each log(EsFR) bin for galaxy NGC 628. Black solid lines illustrate the best fitted 
erfc function to the distributions. Right column: Normalized distribution of log(En 2 ) including both detections and non-detections. Non-detections have been 
randomly distributed using a normal distribution with mean and dispersion parameters provided by the fitted erfc function. To help comparison, Gaussian 
functions with the best fitted mean and dispersion are shown by black solid lines. 


and SF surface density is affected by the diffuse emission 
in the tracers used to estimate the star formation (|Liu et al. 
2011 |Rahman et al.|20lT]|Calzetti et al.|2012[|Momose et al' 

2013 i. In their study of the spatially resolved star formation 
law in the galaxies NGC 3521 and M51, Liu et al. ( |201 1| ) sug¬ 
gest that the difference observed in the exponents of the E H 2 - 
Esfr power law relation in different works (namely, linear vs. 
super-linear relation) is the result of removing or including 
the local diffuse background. 

In order to investigate the effect of the fraction of Esfr not 


associated with current star formation on our Bayesian fitting 
technique, we apply a simple approach to subtract a cirrus 
component of the 24/rm map. A detailed description of our 
methodology is presented in Appendix [A| Our approach con¬ 
sists basically of generating a smoothed map of the 24/mi 
image using a gaussian kernel, which is subsequently sub¬ 
tracted from the original 24//m map. The background sub¬ 
tracted 24/rm map is then used to estimate the Esfr using 
Equation [4] 

With the 24/rm background subtracted Esfr measurements 
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Log(E H2 (M sol pc 2 )) 
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Fig. 17.— Effect of including non-detections in the Bayesian fitting of £h 2 vs. £sfr relations for the grid sample. As in Figure |~j~5] NGC 6946 is in the left 
panel, NGC 628 in the middle and M101 in the right. Black dots represents only detections, i. e., grid points with £h 2 > 10 Mq pc -2 . Red dots show non-detection 
grid points. Solid lines are built by selecting the peak values of each parameter distribution of the relation log( M ^-f R kpc - 2 ) = A + a log( 5Q ^ m pc - 2 ) + Cscat (see 
Figure [19) . The black solid line illustrate the fit to detections only, the red solid line shows the fit to the grid points with non-detections included. 


« 

CT* 

o 


- Fit £ h2 detections only 

- Fit E h2 with non-defections 


NGC.,6946 


-4 


Fig. 18. 


- Fit S H2 detections only 


NGC.'6 2 8 



- Fit S H2 detections only ^ -j Q y' 

- Fit £ h2 with non—defections y 


0 1 2 
Log(S H2 (M so | pc" 2 )) 


0 1 2 
Log(E H 2(M s<> | pc" 2 )) 


0 1 2 
Log(E H2 (M sol pc" 2 )) 


Same as Figure fTT] but Ssfr measurements have been corrected by a cirrus component in the 24//m map. 

in hand, we perform the same a pproac h to distribute the non¬ 
detections presented in Section 4.3.3| The resulting Eh 2 vs. 


Esfr scatter plot for the three galaxies including both detec¬ 
tions and non-detections are shown in Figure |T8] Again, we 
applied our Bayesian analysis of the K-S relationship to the 
case considering E H 2 detections only and to the case when we 
include both detections and non-detections. Figure [l9| shows 
the probability distributions of the parameters a, A ana cr from 
the Bayesian regression fit for each case and for each galaxy 
in our sample. Table [7] shows the peak and 90 % HDI of the 
probability distributions for each fitting parameter. The rela¬ 
tions for each galaxy based on using the peak values of each 
parameter distributions are shown as lines in Figure p~8] 

We do not find a significant change in the distribution of the 
slope between the cases when a cirrus component of the 24//m 
image is subtracted from the Esfr measurements (th is Sec - 
tion) and when no correction has been made (Section |4.3.3[ i. 
Although the slope distributions are slightly shifted to larger 
values when the cirrus correction is applied, the overlapping 
region of the distributions is significant between the two cases. 
On the other hand, the coefficient A is shifted to smaller val- 


but it can move the coefficient A to smaller values due to the 
migration of grid elements to lower values of Esfr. 

5. DISCUSSION 

5.1. Molecular cloud properties 

Based on the broad parameter distributions generated by 
the Bayesian regression fit shown in Figure [8] we do not find 
strong evidence for a R - cr,. relation for the CO( 1-4 0) com¬ 
plexes in the three galaxies studied here. Higher resolution 
observations do not show evidence for a R-o v relation either, 
as the same conclusion can be derived for CO(2 —4 1) clouds 


from Figure 
galaxies M5 


m I, 

IJVO 


In their comparative study of GMCs in the 

_ _ 33 and the Large Magellanic Cloud,[ Hugh es] 

et ah (|2013ji arrive at the same conclusion. They find that, 
when the CO emission is decomposed into GMC-like struc¬ 
tures, the correlation between size and line width is weak or 
almost nonexistent in individual galaxies, or even in the ag¬ 
gregate sample smoothed to a common resolution. 

Studies of the properties of resolved structures have pro¬ 
vided some insights on the nature of the turbulence in the 
ISM and its role in the star formation process both in the 


ues when the 24 pm cirrus component is subtracted from the 

Milky Way (Heyer et al. 2009. Solomon et al. 

vO 

OO 

Shetty 

Esfr measurements for the galaxies NGC 6946 and NGC 628. 

et al. 

2012 

and in nearby galaxies (Bolatto et a 

2008 

Sheth 

Thus, the subtraction of the cirrus component of the Esfr val- 

et al. 

2008 

Rebolledo et al. 2012). Several simu 

lation 

works 


ues has little or no impact on the distribution of the slopes, have proposed a scenario in which the turbulence observed in 
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TABLE 7 

Bayesian regression parameters for £ H 2 -Ssfr relation using an hexagonal grid. 


Grid elements considered in the fit 

Q 

90% HDI 

A 

90% HDI 

(7 

90% HDI 

NGC 6946 

Eh 2 detections only 

0.63 

[0.48, 0.68] 

-1.57 

[-1.61,-1.54] 

0.22 

[0.19, 0.23] 

Eh 2 non-detections included 

0.38 

[0.28,0.41] 

-1.59 

[-1.63,-1.56] 

0.23 

[0.21,0.24] 

Eh 2 detections only 24/j,m background subtracted 

0.68 

[0.53, 0.78] 

-1.63 

[-1.68,-1.60] 

0.25 

[0.21,0.26] 

Eh 2 non-detections included 24/im background subtracted 

0.38 

[0.26,0.41] 

-1.65 

[-1.70,-1.62] 

0.26 

[0.24, 0.27] 

On-arm sample 

0.63 

[0.50, 0.74] 

-1.45 

[-1.50,-1.43] 

0.21 

[0.18,0.22] 

Inter-arm sample 

0.23 

[0.08, 0.33] 

-1.75 

[-1.81, -1.73] 

0.10 

[0.06, 0.11] 

On-arm sample 24[im background subtracted 

0.73 

[0.55, 0.82] 

-1.51 

[-1.56,-1.48] 

0.23 

[0.19, 0.24] 

Inter-arm sample 24 (im background subtracted 

0.28 

[0.06, 0.42] 

-1.87 

[-1.94,-1.84] 

0.12 

[0.08,0.15] 

NGC 628 

£ H 2 detections only 

0.63 

[0.44, 0.74] 

-1.87 

[-1.92,-1.83] 

0.23 

[0.20, 0.24] 

Eh 2 non-detections included 

0.63 

[0.52, 0.65] 

-1.85 

[-1.91,-1.83] 

0.20 

[0.18,0.21] 

Eh 2 detections only 24/im background subtracted 

0.73 

[0.53,0.91] 

-1.97 

[-2.03,-1.91] 

0.28 

[0.25, 0.30] 

Eh 2 non-detections included 24/j,m background subtracted 

0.68 

[0.58, 0.75] 

-1.97 

[-2.02,-1.93] 

0.26 

[0.23, 0.27] 

On-arm sample 

0.68 

[0.43, 0.85] 

-1.79 

[-1.86,-1.74] 

0.24 

[0.20, 0.26] 

Inter-arm sample 

0.38 

[0.16,0.52] 

-2.03 

[-2.10,-1.97] 

0.17 

[0.14, 0.19] 

On-arm sample 24fim background subtracted 

0.78 

[0.49, 1.00] 

-1.89 

[-1.97,-1.82] 

0.28 

[0.24, 0.31] 

Inter-arm sample 24(im background subtracted 

0.48 

[0.19,0.71] 

-2.15 

[-2.23, -2.07] 

0.22 

[0.18,0.25] 

M101 

£ H 2 detections only 

0.98 

[0.64, 1.23] 

-1.69 

[-1.79,-1.59] 

0.44 

[0.40, 0.47] 

Eh 2 non-detections included 

0.83 

[0.69, 0.86] 

-1.67 

[-1.73,-1.61] 

0.30 

[0.27, 0.30] 

Eh 2 detections only 24/zm background subtracted 

1.03 

[0.70, 1.33] 

-1.71 

[-1.82,-1.61] 

0.47 

[0.42, 0.50] 

£h 2 non-detections included 24^,m background subtracted 

0.83 

[0.69, 0.88] 

-1.71 

[-1.78,-1.65] 

0.33 

[0.30, 0.33] 

On-arm sample 

0.88 

[0.46, 1.20] 

-1.53 

[-1.65,-1.41] 

0.48 

[0.43, 0.54] 

Inter-arm sample 

0.48 

[0.18,0.68] 

-2.09 

[-2.18, -2.00] 

0.19 

[0.16, 0.21] 

On-arm sample 24^m background subtracted 

0.93 

[0.50, 1.27] 

-1.55 

[-1.68,-1.42] 

0.51 

[0.45, 0.56] 

Inter-arm sample 24(im background subtracted 

0.53 

[0.23, 0.78] 

-2.11 

[ -2.23, -2.04] 

0.21 

[0.17,0.23] 


molecular clouds originates on the largest scales in the ISM 
( Heyer & Brunt|2004J |Brunt et al.|2009) . According to this 
framework, the velocity dispersion will remain constant for 
scales larger than the size scale driving the turbulence in the 
ISM ( Ossenkopf & Mac Low| 2002j >. In their study of the 
dens e T5M in the Central M olecular Zone (CMZ) of the Milky 
Way, She tty et al.| ( |2012j > find that the trend observed in the 
cr v -R relation is independent of the presence of dense struc¬ 
tures in the region probed by the observations. They interpret 
this result as evidence in favor of the turbulence being injected 
at larger scales in the CMZ. Over the size scales probed by 
their observations (R ~ 2-40 pc), they do not find a clear 
flattening of the a V —R relation, leading them to suggest that 
the turbulent velocities in the CMZ are driven on scales larger 
than 30 pc. 

The <7,. — R scatter plots of the identified structures in the 
three galaxies in our sample may provide some insights on 
the scales at which turbulence is driven in the ISM in spi¬ 
ral galaxies. Among the three galaxies in our sample, NGC 
628 presents the least inclined disk, offering the best place 
to estimate the cloud velocity dispersions without the influ¬ 
ence of cloud blending or crowding. The a v -R scatter plots 
for complexes (Figure [8]> and for clouds (Figure |9]> show 
that, for the scales probed by our observations in NGC 628 
(~ 70-200 pc), the velocity dispersions reach a maximum 
value ~ 7 km s _l . A similar trend is found for complexes and 
clouds for M101. However, given the large uncertainties in 
the slope due to the limited dynamic range in size, as shown 
in Tables [4] and [5] it is difficult to draw any conclusion from 
the observed a v —R relation about the nature of the turbulence 
in the ISM in NGC 628 and M101. 

In the case of NGC 6946, we observe some complexes and 
clouds with velocity dispersions ~ 10 kms -1 . In Paper I, 


our observations of the eastern part of the disk of NGC 6946 
showed that the regions with the most active star formation 
are the regions with the largest clouds and complexes, and 
with the largest velocity dispersions as well. Moreover, in 
one of the active star forming regions of NGC 6946, we find a 
large velocity gradient in the CO emission map. Thus, the dif¬ 
ference in the star formation rate in those regions, along with 
other kinematic features in the gas, may be influencing the re¬ 
lation between <j v and size with respect the other regions of 
the disk. In fact, in the center of NGC 6946|Donovan Meyer| 
et al.| (|2012|) find a set of GMCs with velocity dispersions 
larger than other clouds with similar sizes. These clouds may 
be dynamically affected by the presence of a nuclear bar in 
NGC 6946, resulting in larger line widths compared to com¬ 
plexes or clouds with similar sizes located in other regions 
of the disk. A complete census of resolved GMCs in several 
environments such as the center of a barred galaxy, regions 
located in a strong spiral arm, and regions with high and low 
star formation activities is needed to draw a robust picture of 
the influence of those environments on the turbulent nature of 
the molecular clouds. 


5.2. Kennicutt-Schmidl relation 


As 


et al. 


has been discussed recently (|Liu et al.||2011| |Calzetti 
2012| Momose et al. 201 3|>, deriving the tneK-S relation 


and its intrinsic uncertainty is a non-trivial task. It depends 
on several factors such as the selection of the tracers used to 
estimate the gas content and the SFR activity over a particular 
region, the statistical tools to find the best fit to the power 
law, or the physical scale probed by the observations, among 
others. Analyses that combine different approaches to deal 
with the effects of these factors may drive the estimated fit 
parameters and uncertainties to a wide range of values. 

While our analysis is subject to some of the same limi- 
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Fig. 19.— Probability distributions of the parameters a, A and a from the Bayesian regression fit to the relation log(- 


^SFR 


— )=A+olog(- 


£h2 


r? )+6scat- 


v Mq yr 1 kpc 2 2 ° 50 Mq pc 2 J 

Top row panels show the parameter distributions for NGC 6946. middle row panels show the distributions for NGC 628, and bottom row panels show the 
distributions for M101. Left column panels show the distributions of the slope a, middle column panels show the distributions of the intercept coefficient A, and 
the right column panels show the distribution of the dispersion of the scatter, a. Black histograms show the resulting distributions considering only S R 2 detections 
in the Bayesian regression fit. Blue histograms illustrate the resulting distributions with non-detections included. Light blue histograms show the distributions 
considering only detections, but with Esfr values corrected by the 24 /a m cirrus component. Finally, red histograms show the distributions considering points 
with Sjj 2 non-detections included and Esfr values corrected by 24//m cirrus component. 


tations, our Bayesian approach provides us with more real¬ 
istic estimates for uncertainties of the regression fit param¬ 
eters. Among our sample of galaxies, NGC 6946 shows 
the most constrained K-S relationship at the scale probed by 
CO(l -A 0) complex boundaries. On the other hand, the pa¬ 
rameter distributions for NGC 628 and M101 are poorly con¬ 
strained, particularly in the distribution of the slope. The lim¬ 
ited dynamic range present in the CO observations of NGC 
628 on the one hand, and the large variation of the SFR ob¬ 
served in M101 on the other, makes the derivation of a power 
law unstable for these two galaxies. The increased resolution 
offered by the CO(2 —>• 1) clouds is translated into parameter 
distributions consistent with super-linear relations for NGC 
6946 and M101, but with increased uncertainties. Addition¬ 
ally, the dispersion of the intrinsic scatter, a, shows a factor 
of ~ 2 higher values compared to the CO( 1 —► 0) complexes 
for the three galaxies. This is consistent with the scenario pro¬ 
posed by Calzetti et al.|( |2012) , where the slope and the scatter 
of simulated K-S relations depend on the size scale probed by 
the region considered in the analysis. 

When we study the K-S relation using hexagonal grids, on 
the other hand, we find slope distributions consistent with 


sub-linear relations for NGC 6946 and NGC 628, and a less 
constrained slope distribution for M101 that covers sub-linear 
and super-linear relations. Given that the grid approach al¬ 
lows us to include Eh 2 non-detections in the K-S analysis, we 
are able to assess the change in the fitting parameters when 
non-detections are considered. Once the non-detections are 
included in the fit, we find that the slope distributions for all 
three galaxies are consistent with sub-linear relations. 

Including Ejn non-detections can clearly have an impact on 
the power law index of the K-S relationship. The left column 
of Figure [19] show that in the case of NGC 6946, the distri¬ 
bution of the slope a is shifted to smaller values when E H 2 
non-detections are included. On the other hand, for M101, 
although the slope distribution of the K-S relationship consid¬ 
ering Eh 2 detections and non-detections is shifted to values 
smaller than the case with detections only, both distributions 
significantly overlap. For NGC 628, both slope distributions 
are consistent. The overall distribution on the Esfr axis of 
the E H 2 non-detections ultimately drives the effect of includ¬ 
ing them in the linear fit. In the case of NGC 6946, ~ 75% 
of the E H 2 non-detections have E SFR > 0.006 M 0 yr -1 kpc -2 . 
On the other hand, this percentage decreases to 34% in the 
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Fig. 20.— Radial distributions of Eh2 (left), Esfr (center) and Esfe (right) of the regions observed in NGC 628, including both CO(l —> 0) complexes 
and CO(2 > 1) clouds. Red dots represent complexes, while black dots illustrate the clouds. While the molecular gas surface density Eh 2 shows a nearly flat 
distribution over the disk, the SFR surface density Esfr expands a range of several order of magnitudes. Thus, the star formation efficiency Esfe shows a range 
of values from 5 X 10 5 to 0.006 Myr 1 , but with no clear peaks as we observed in NGC 6946 in Paper I. 





Distance to galactic center (kpc) Distance to galactic center (kpc) Distance to galactic center (kpc) 


Fig. 21.— Radial distributions of Efe> (left), Esfr (center) and Esfe (fight) of the regions observed over the southern arm of M101, including both CO(l > 0) 
complexes and CO(2 > 1) clouds. Red dots represent complexes, while black dots illustrate the clouds. Similarly to what we have found in galaxy NGC 6946, 
clear peaks of Esfr and Esfe are observed in M 101. 


case of NGC 628, and to 29% in the case of M101. Thus, a 
large fraction of grid points with undetected CO but high SFR 
bias the slope of the linear fit to smaller values in the case of 
NGC 6946, while a large fraction of E H 2 non-detections with 
low SFR keeps the fit parameters similar to the trend found 
when only Em detections are considered in the case of NGC 
628. 

The larger fraction of grid elements with high SFR but non- 
detected CO in NGC 6946 compared to NGC 628 and M101 
could be related to several observational effects. For exam¬ 
ple, spatial offsets between CO emission and regions of active 
SFR may be due to cloud evolution. Thus, maps with enough 
angular resolution to resolve this spatial offset will be able to 
populate the low-CO and high-SFR region of the K-S rela¬ 
tion. In our grid analysis, we have used a fixed grid size of 8" 
for the three galaxies, which translates into different physical 
scales due to differences in the galaxy distances. For NGC 
6946, the grid size is 213 pc, while for NGC 628 and M101 
the grid size is 283 pc and 287 pc respectively. Given the 
similar grid sizes of the three galaxies, we conclude that, if 
cloud evolution were responsible for the grid elements with 
high SFR but weak CO in NGC 6946, similar trends would 
have been detected in NGC 628 and M101 as well. 

In contrast to recent studies, we do not observe a significant 
change in the slope distribution when a 24//m cirrus compo¬ 


nent is subtracted from the Esfr measurements. In Figure 18 


we notice that this correction has migrated some of the Esfr 
values to the region in log(EsFR) < -2.5, especially for NGC 
6946. For NGC 6946 and NGC 628, the net effect of the 
24//m cirrus component subtraction is to move the intercept 
coefficient A to smaller values with respect to the case where 
no correction is applied. For M101, no detectable change has 
been observed. 

Alternatively, regions of low CO brightness and high SFR 
may be related to variations of the Xqo factor. By assuming 
a constant value for .Xco across the surveyed area, we may 
be underestimating the amount of molecular mass over re¬ 
gions where CO is dissociated by the radiation field, but FF 
is still present. Direct comparison with similar spatial resolu¬ 
tion maps of other mass tracers such as infrared dust emission 
will provide independent constraints on the amount of total 
molecular mass in these regions. 

5.3. Distribution of the star formation properties in the disk 

In Figures[20]andl2T]we show the distribution of Eh 2 , Esfr 
and Esfe (oFE is tKe star formation efficiency defined as 
Esfe = Esfr/Sh 2 ) with respect to the galactocentric radius 
for the complexes and clouds identified in galaxies NGC 628 
and M101 respectively. We refer the reader to Paper I (see 
Figure 14) for NGC 6946 plots. In the case of NGC 628, 
we observe a roughly uniform distribution of molecular gas 
density, with no clear trend with radius. There are a couple 
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Fig. 22.— Radial distributions of the Fourier amplitude modes I m normalized to the Iq mode for galaxies NGC 6946 (left), NGC 628 (center) and M101 
(right), using the 3.6/xm images from SINGS. Vertical dashed lines illustrate the galactocentric radial extent covered by our CO(l —¥ 0) observations. The radial 
distribution of the angular phase for m = 2 (02) is shown in the lower panels. NGC 628 shows prominent I 2 mode in the region we observed the molecular gas, 
and 02 distribution is consistent with a two-arm spiral structure of the disk. On the other hand, higher Fourier modes are also important for the multi-arm galaxies 
NGC 6946 and M101. Although I 2 is important in the center of NGC 6946, the flat distribution of 0 2 reflects the presence of a bar structure, which is different 
from a two-arm structure as seen in NGC 628. 


of CO(2 —» 1) clouds with Eh 2 >100 M 0 pc~ 2 which are lo¬ 
cated on on-arm regions. The star formation surface density, 
on the other hand, presents a wider distribution of values rang¬ 
ing from ~ 1 0~ ' to ~ 5 x 10 _1 M 0 yr 1 kpc 2 , resulting in a 
wide distribution of SFE surface densities as well. 

In the case of M101, the molecular gas surface density dis¬ 
tribution is similar to that observed in NGC 628, i. e., there 
is not a clear trend of Eh 2 with galactocentric radius. On the 
other hand, the Esfr shows pronounced peaks at galactocen¬ 
tric radius ~ 7, 10 and 17 kpc, which translate to peaks in 
the SFE for those regions as well. Those regions have been 
highlighted in Figure [2] A similar behavior of the Esfr was 
observed in NGC 6946 in Paper I, where we identified two 
specific regions in the eastern part of the disk with high SFE. 

Morphologically speaking, NGC 6946 and M101 seem to 
be different from NGC 628. Whereas NGC 6946 and M101 
show a multi-arm spiral structure, NGC 628 shows a two-arm 
spiral pattern in the inner part of the disk (Foyle et al .|2010| >. 
In order to quantify the significance of grand design spiral 
structure in the three galaxies in our sample, we use the same 
approach that we applied in Paper I for the galaxy NGC 6946 
to define the spiral arms. We estimate the arm amplitude by 
performing a Fourier series expansion of the old stellar popu¬ 
lation brightness distribution traced by the 3.6 (i m map from 
SINGS. Firstly, for each galaxy we have deprojected the 3.6 
fim images using the parameters shown in Table [Tj Then, we 
decompose the images into radial bins with a width of 3"75, 
overlapped by 50% of the width for smoothness of the decom¬ 
position. For every radial annulus, we fit the function given 
by 


f^max 

I(r,4>) = / 0 (r) + y^/ m (r)cos[m x (</>-</>„,(r))] (14) 

m-l 

where /„, is the arm amplitude and </>,„ is the phase for each 
Fourier component m , and rn max is the maximum component 
considered in the expansion. We have used the MPFIT pack¬ 
age (Mar kwardt|2009] > to determine the best-fit parameters in 
Equation ( [T4| . The relative importance of the individual com¬ 
ponents is ultimately determined by the complexity of the arm 
structure in the galaxy. In the case of grand-design galaxies, 
the spiral arm structure is well recovered by using components 
up to m = 4, with a predominance of the mode m = 2 (|Foyle| 


et al.|2010). Flocculent and multiarm galaxies, however, usu¬ 
ally require higher components to recover their complex arm 
structure. In Figure [22] we show the resulting /,„ amplitudes, 
normalized to the amplitude of m = 0, of the three galaxies 
as a function of the galactocentric radius. Additionally, Fig- 
ure|22|shows the phase for m = 2, (tMr), for the three galaxies. 
NGC628 is the galaxy that shows a distribution of amplitudes 
that resembles the amplitude spectrum expected for a two-arm 
grand design galaxy, i. e, a prominent mode m = 2 in most of 
the disk traced by 3.6/rm, and <bi(r) steadily increasing with 
radius. On the other hand, for galaxies NGC 6946 and M101 
the m = 2 amplitude is not the predominant mode all across the 
disk, but higher modes are also significant. Although NGC 
6946 presents a prominent h in the center of the disk, 02 is 
approximately constant with radius in the inner part. This be¬ 
haviour of fair) is consistent with the presence of a bar in the 
center of NGC 6946, and morphologically different from the 
two-arm spiral structure observed in NGC 628. 

Why does the distribution of star formation efficiency seem 
to be different for some regions of flocculent/multi-arm spiral 
galaxies compared to what is observed in spiral galaxies with 
symmetric spiral arm structure? Several authors have used 
hydrodynamic simulations to model flocculent and multi-arm 
spiral galaxies by using a “live” stellar disk potential (Dobbs 
& Bonnell|2008[ |Wada et al.|2011[ ). They find that the struc- 


tures observed in these types of galaxies are driven by mecha¬ 
nisms different from the stellar density waves present in grand 
design spiral galaxies. In their simulations, the gas flows into 
the local potential minima, where it stays and condenses to 
form the dense gas that triggers the star formation. Thus, the 
gas traces the potential minimum, and no clear offset between 
the gas tracers and spiral potential is observed in the simu¬ 
lations. Nevertheless, maybe the most remarkable aspect of 
these simulations is that the spiral arms are not steady on time 
scales of a few 100 Myr (Clarke & Gittins] [2006) >, with fre¬ 
quent collision and merging of the spiral arms. Moreover, the 
simulations show that the largest GMCs are located in the re¬ 
gions where the spiral arms collide. The collision of spiral 
arms may induce the formation of gas overdensities, trigger¬ 
ing star formation. In our observations, both NGC 6946 and 
M101 show spot regions of high star formation on the spi¬ 
ral arm (Figure [2T] and Figure 14 in Paper I), which may be 
regions where previous spiral arms have collided. 
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Fig. 23.— £h 2 vs - ^sfr relations using a uniform grid, but in this case we have included the classification of on-arm and inter-arm regions yielded by the 
Fourier decomposition method. Black dots represents the grid elements located in on-arm regions. The red dots correspond to inter-arm regions. Black solid lines 
illustrate the regression fit to the on-arm grid elements, while the red solid lines are the regression fit lines for the inter-arm regions. Black dotted lines represent 
the regression fit performed to the full sample of grid elements. We include the histograms of the SFR surface density on the right panels, along with the number 
of grid elements classified as on-arm or inter-arm regions. 


For NGC 628, on the other hand, the maximum values of 
SFE are located across a range of radii in the disk (Figure |20|i. 
This is consistent with the picture of the formation of massive 
molecular clouds attributed to the passing of spiral density 
waves in gran d-design galaxies ( |Kim & Os triker 20021 Wada 
|& K oda 2004). In this scenario, regions of high star formation 
activity will be more continuously distributed in radius as they 
are aligned with the spiral structure, in contrast to the situation 
observed in NGC 6946 and M101. 

5.4. On-arm us. inter-arm £h2-5jsfr re ^ at ‘ on - 

The Fourier decomposition analysis of the 3.6//m image al¬ 
lows us to identify the regions associated with the arms and 
those associated with the inter-arm zone. We have classified 
the on-arm and inter-arm regions following the same approach 
adopted in Paper I. A reconstructed image of the 3.6//m is 
generated from all the modes considered in the Fourier de¬ 
composition. For each radial bin, we have selected the 35% 
highest pixels to generate the mask representing the on-arm 
regions. That mask is then applied to the grid map in or¬ 
der to assign grid elements to on-arm or inter-arm regions. 
For simplicity, we have only considered £ H 2 detections in 
our on-arm vs. inter-arm analysis. In Figure [23] we show 
the £h 2 -Ssfr relations for the three galaxies in our sample 
including the on-arm vs. inter-arm classification. For NGC 
6946 and M101, it is clear that on-arm regions are distinct 
from inter-arm regions, consistent with what we found in Pa¬ 
per I for complexes and clouds. The on-arm regions present 
higher SFR than inter-arm regions for similar molecular gas 
surface densities. In NGC 628, although the main group of 
points contains both on-arm and inter-arm clouds, we observe 
a group of on-arm grid elements that deviates toward higher 
SFR surface density. 

In order to investigate whether the K-S relation changes re¬ 
spect to the location in the disk, we have applied the Bayesian 
regression fitting approach (Equation ED to the on-arm and 
inter-arm regions. Figure [23] shows the resulting fitting lines 
for the three galaxies in our sample built from the peak value 
of each parameter distributions. In Table [7] we present the 
peaks and 90 % HDI values of the parameter distributions 
from the Bayesian regression fit. The overall result obtained 
for the three galaxies in our sample is roughly similar: the dis¬ 
tributions of slopes of inter-arm regions are shifted to smaller 
values than those found in K-S relations for on-arm regions. 


Also, the distributions of slopes considering on-arm regions 
are similar to those found when we consider the full sample of 
grid elements. Regarding the distribution of the intercept co¬ 
efficient A, we observe that this parameter reaches the smaller 
values for inter-arm grid regions, followed by the case when 
all the grid points are considered. We conclude that low level 
SFR regions mainly present in inter-arm regions does not af¬ 
fect significantly the slope of the K-S relation when all the 
grid values are considered. In contrast, the K-S linear rela¬ 
tions for the full sample of grid points are biased to smaller 
values of SFR respect to the on-arm K-S relations, which is 
evident from the smaller values of the A parameter yielded by 
the Bayesian regression fit for the full sample respect to the 
on-arm regions. 

The different behaviours found for on-arm vs. inter-arm 
regions presented above may be reflecting some underlying 
physical process in the formation of stars in spiral galaxies. 
For instance, in their study of the barred spiral galaxy NGC 
4303, |Momose et ak| ( ]2010| found differences in the SFE be¬ 
tween the spiral arms, the bar, and the inter-arm regions. The 
fact that inter-arm molecular gas has similar surface density 
compared to on-arm molecular material but with lower SFR, 
may be an indication of the presence of a higher amount of 
denser gas inside the molecular clouds located in the spiral 
arms, especially in the case of M101 where this difference in 
SFR is more prominent. The spiral shocks (maybe present on 
NGC 628), or convergence of gas flows (maybe more related 
to galaxies NGC 6946 and M101) may lead to an increase in 
the amount of dense gas in spiral arms. Several authors have 
investigated the nature of the extragalactic scaling relation of 
SFR and dense molecular gas (|Gao & Solomon 2004] |Lada| 
■12012)) . |Lada et ah (2012 1 suggest that total star torma- 
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tion rate in a molecular cloud is linearly proportional to the 
amount of dense gas present in the cloud. Thus, the steeper 
£ H 2 -£.sfr relation observed in Figure [23] for the on-arm re¬ 
gions, more prominent in NGC 6946, may indicate the pres¬ 
ence of denser gas created by processes that depends on the 
global properties of the dynamics of the gas. 

In our analysis using the uniform grid, we are including 
regions from the inter-arm zone where the contribution from 
the diffuse component unrelated to the current star formation 
may bias the fitted linear relation to smaller slopes. The effect 
of the diffuse component is to overestimate the SFR in the 
lower end of the distribution, flattening the slope of the E H2 - 
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Same as Figure |23] but a subtraction of a 24/un smooth component is applied to Esfr. 

Esfr linear relation in the logarithmic space. This may be 
the cause of the flatter relation we find in our uni form g rid 
analysis compared to the relations found in Section |4.3. 1 |for 
complexes and clouds. As we defined our regions based on 
CO(l —> 0) complexes and CO(2 —»• 1) clouds, we are likely 
selecting regions less affected by the diffuse contribution. The 
small impact of the diffuse emission in the K-S law in regions 
of high surface density was already reported byjRahman et al.| 

(120111) in their study on the central region of the disk of NGC 

425 ^ 


Following the approach already introduced in Section 4.3.4 
and detailed in Appendix [A] we have subtracted a diffuse 
component in the 24//m not associated with current formation 
from the Esfr measurements. Figure 24 shows the result¬ 


ing K-S relations for the three galaxies in our sample, with 
the Esfr measurements corrected by the cirrus component 
in the 24//m. We notice that differences between on-arm vs. 
inter-arm regions is similar to the case when no cirrus correc¬ 
tion is applied, with on-arm regions presenting higher Esfr 
than inter-arm regions, being more evident for NGC 6946 and 
M101. Therefore, our approach to subtract a diffuse compo¬ 
nent in the 24//m map has not introduced a significant change 
in the on-arm vs. inter-arm region differentiation. 

We apply the Bayesian regression fitting to the K-S relation, 
but now with the Esfr measurements corrected by the cirrus 
component in the 24 pm. Table [7] shows the resulting peaks 
and 90 % HDI values of the parameter distributions. The over¬ 
all distribution of the slopes is consistent with the case where 
no subtraction of the 24 pm smooth component is applied. In 
the same direction, the distributions of the intercept coeffi¬ 
cient A and the dispersion of the intrinsic scatter er are con¬ 
sistent with the case when Esfr is not corrected for the 24//m 
diffuse component. In the case of NGC 6946, the Bayesian 
regression fitting to the K-S relation for on-arm regions yields 
larger slopes than inter-arm regions, which is consistent to the 
case with no cirrus correction. On the other hand, for NGC 
628 and M101 the difference between on-arm vs. inter-arm 
region is less prominent. Again, this is consistent to the K-S 
relation differences between the on-arm vs. inter-arm regions 
when no correction of the 24/im diffuse component is applied 
to the the Esfr values. 

Although our analysis of on-arm and inter-arm regions only 
includes Eh 2 detections, including the Eh 2 non- detecti ons 
would not change our main result. In Section |4.3.3| we 
see that, in general, the K-S relations present flatter slopes 
when we include Eh 2 non-detections. Because the Eh 2 non¬ 
detections have a stronger effect on inter-arm grid points, our 


finding that inter-arm regions have K-S relations with flatter 
slopes than on-arm regions would be preserved. 


5.5. Caveats and uncertainties of the present work 

The coverage of our observations have been limited to re¬ 
duced regions of the molecular gas disk of NGC 6946, NGC 
628 and M101. Therefore, our conclusions can only be ap¬ 
plied to the regions covered by our observations and cannot 
be interpreted as global properties of the galaxies included in 
our sample. A full census of the GMCs population across the 
disks of nearby galaxies is needed in order to perform a de¬ 
tailed study of molecular cloud properties and their relation 
with the surrounding environment. Nevertheless, high sensi¬ 
tivity maps with resolutions ~ 60 pc that cover a significant 
part of the molecular disk are observationally challenging, al¬ 
though recent studies have achieved this goal in M51 ( Schin-| 
nerer et al.|2013| l. In that direction, we have also successfully 
completed observations of the full disk of NGC 6946 with 
CARMA (Rebolledo et ah, in preparation). With the help of 
facilities such as ALMA, we will be able to obtain high qual¬ 
ity images at resolutions close to GMC sizes, allowing us to 
perform unbiased analyses of the physical properties of the 
molecular gas in nearby spiral galaxies. 

Equations |4]and[5]assume a continuous star formation pro¬ 
cess. These prescriptions to estimate the SFR may break 
down for scales close to GMC sizes, as we may be isolating 
regions with a single stellar population at a specific age, espe¬ 
cially with our CG(2 —> 1) observations. In order to investi¬ 
gate the age sensitivity and the intrinsic scatter in the SFR es¬ 
timates by using Equations[4]and[5]at small scales, [Leroy et al.| 
( |2012j ) simulate the Ha and FUV emission after an instanta¬ 
neous burst of star formation. They show that Ha is mostly 
emitted in the first 10 Myr. On the other hand, FUV presents 
significant emission up to 65 Myr after the burst. This model 
predicts a factor of ~2 uncertainty inferring the SFR from Ha, 
and a factor ~3-4 uncertainty inferring SFR from FUV. 

Additionally, these equations assume that the ionizing star 
clusters are massive enough in order to have the high mass tail 
of the IMF well populated (|Kennicutt||1998| |Kennicutt et al.| 
2007|l . However, this is likely not to be true in regions of 
low Esfr, where the SFR tracer emission could be dominated 
by low mass stars. In an effort to understand the different 
sources of the observed scatter in the K-S relation in spatially 
resolved studies,[Kruijssen & Longmore|(j2014| present a for¬ 
malism based on a simple uncertainty principle for SF. They 
show that for an idealized galactic environment typical of a 
disk galaxy, the minimum size required to recover an under- 
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lying SF relation is ~ 500 pc, and is determined by the ability 
to sample independent star-forming regions. This is a factor of 
~3-5 larger than the sizes sampled by our CO( 1 —► 0) obser¬ 
vations, a factor ~6-10 larger than the CO(2 —> 1) cloud sizes, 
and ~2 larger than the area of the grid elements. Thus, the 
high scatter observed in our derived K-S relations, especially 
for CO(2 —> 1) clouds, may be the result of our incomplete 
sampling of independent star forming regions over the areas 
covered by our observations. Alternative approaches to mea¬ 
sure the SFR at small scales (for instance, via observations of 
resolved star clusters) can be compared to the SFR traced by 
FUV or Ha, allowing us to assess the level of uncertainties in 
our methodology to estimate Esfr- 

6. SUMMARY 

We have performed high resolution observations of the 
molecular gas in the nearby spiral galaxies NGC 628 and 
M101. Using CARMA, we have observed the CO(l —> 0) 
emission line over regions of the disk with active star forma¬ 
tion and offset from the galactic center. These observations 
have supplemented the CO( 1 —>• 0) observations of the north¬ 
eastern region of NGC 6946 reported in Paper I. Higher reso¬ 
lution observations of CO(2 —► 1) toward the brightest regions 
observed in CO(l —» 0) have allowed us to resolve some of the 
largest GMCs. The results are summarized as follows: 

1. Usi ng the cloud-finding algorithm CPROPS (Section 
|4. 1[ ), we have identified 112 CO emitting complexes 
with typical sizes of ~ 150 pc. The higher resolution ob¬ 
servations of CO(2 —► 1) towards the brightest regions 
detected in CO(l —► 0) allowed us to find 144 structures 
with sizes ~ 70 pc which we have identified as GMCs. 
Properties such as size, line width and luminosity were 
derived using the moment approach implemented in the 
CPROPS package. 

2. We observed that the size-line width relations for the 
CO(l —> 0) complexes and CO(2 —> 1) clouds present 
significant scatter, with some of the regions in NGC 
6946 presenting excess velocity dispersions. On the 
other hand, line widths in NGC 628 are limited to 
< 7 km s -1 . Additionally, our M vlr -Lco scatter plot 
for both complexes and clouds are consistent with our 
choice of the CO to H 2 conversion factor Xqq=2 x 
10 20 cm _2 (K km s -1 ) -1 . 

3. Linear fitting to the scaling relations between the re¬ 
solved properties of the identified structures is per¬ 
formed using a Bayesian regression approach. In agree¬ 
ment with recent studies of GMCs in nearby galax¬ 
ies, we do not find strong correlations in the size-line 
width relations when resolution bias correction is ap¬ 
plied. The reduced dynamic range in the sizes of the 
identified structures does not allow to invoke a physical 
explanation to the apparent presence of an upper limit 
in the velocity dispersions observed in NGC 628. 

4. The Bayesian linear regression for the E H 2 - Esfr re¬ 
lation for complexes is consistent with a super-linear 
K-S relation in the case of M101. On the other hand, 
the lower tail of the slope distribution reaches values 
below 1 for NGC 6946 and NGC 628. In general, the 
slopes for CO(2 —y 1) clouds are similar to those found 
for CO(l —> 0) complexes, but with greater uncertainty. 


The distribution of a is consistent with super-linear re¬ 
lations for NGC 6946 and M101, while it reaches val¬ 
ues below 1 in the lower end for NGC 628. The higher 
scatter found in higher resolution maps may be the re¬ 
sult of variations in the star formation activity at local 
scales, which translates to a larger variation of the SFR 
than the molecular gas surface density from region to 
region. 

5. In contrast to the case when we analyse the K-S rela¬ 
tion for structures identified by CPROPS, the Bayesian 
regression fit using a uniform grid yields sub-linear K- 
S relations for NGC 6946 and NGC 628. In the case 
of M101, the slope probability distribution has a peak 
at ~ 1. We do not observe a significant change in the 
slope distributions when a cirrus component in the Ssfr 
is subtracted. When £h 2 non-detections are included 
in the analysis, the slope of the K-S relation may be 
affected depending on the overall distribution of Esfr 
along the y—axis. 

6. Using the 3.6 /jm images of NGC 6946, NGC 628 and 
M101 as stellar mass density tracers, we have imple¬ 
mented a Fourier decomposition approach to identify 
the spiral arm structure present in these galaxies. On- 
arm and inter-arm regions have been defined based on 
the 3.6 /un Fourier reconstructed images. We observe 
that, over the regions observed in the galaxies in our 
sample, the most prominent Fourier mode for NGC 628 
is m = 2, similar to the expected value for a two-arm 
grand design galaxy. On the other hand, NGC 6946 and 
M101 lack a prominent m = 2 mode, showing signifi¬ 
cant higher order modes instead, giving evidence for a 
multi-arm or flocculent structure present in these galax¬ 
ies. 


7. This difference in the distribution of regions with maxi¬ 
mum SFR and SFE may be related to the underlying dy¬ 
namical process that drives the arm structure observed 
in spiral galaxies (Section 5.3 1 . In this scenario, in 
galaxies with nearly symmetric arm shape (e. g., NGC 
628), the spiral shocks are triggering the star forma¬ 
tion along the arms, giving a more uniform distribution 
of SFE peaks across the disk. On other hand, galaxies 
with flocculent or multi-arm spiral structure (e. g., NGC 
6946 and M101) show regions of high star formation 
efficiency on specific regions of the spiral arms (Figure 

f in this paper, and Figure 14 in Paper I), which may 
interpreted as the result of gas flow convergence or 
regions where previous spiral arms may have collided. 
High fidelity observations of the molecular gas and star 
formation covering the full disk of nearby galaxies with 
different morphology will be required to test the sce¬ 
nario proposed here. 


8. The Bayesian regression fit to the Eh 2 _ Ssfr relation 
in logarithmic scale for on-arm and inter-arm regions 
reveals that the distribution of the slopes are shifted 
to higher values for on-arm grid points for the three 
galaxies in our sample. Same picture is obtained when 
the Ssfr values are corrected by the cirrus component. 
Higher star formation activity found in on-arm regions 
may indicate the presence of denser gas, which may be 
induced by processes that depend on the global dynam¬ 
ics of the gas. 
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chosen the size of the gaussian kernel based on the unsharp 
masking approach used by [Rahman et al. (J20 I to estimate 
the diffuse emission on the SF tracer maps in their study of 
the galaxy NGC 4254. They suggest that using kernel sizes 
below ~ 6 kpc overestimates the diffuse component fraction. 
Thus, as a first approximation, we have adopted the 6 kpc as 
the size of gaussian kernel in our approach. This physical size 
corresponds to ~ 170" for NGC 628 and M101, and 225" for 
NGC 6946. The next step is to multiply the smoothed 24//m 
map by the /de factor, and subtract this map from the origi¬ 
nal 24//m map. We select the value of / DE based on the best 
linear relation between the Xsfr from Ha and the cirrus sub¬ 
tracted Xsfr from 24/jm+FUV. We assess the level of linear 
correlation between these two quantities using a Pearson cor¬ 
relation coefficient, r P . We calculate r P for different values 
°f /de = 0.2,0.3,0.4, and 0.6, and we select the value that 
yields the highest r P . Table[8]shows the values for r P for each 
value of / D e- Additionally, we have included the value of r P 
when no smooth background subtraction of 24//m is applied. 
We see that in the case of NGC 6946, /de=0.3 produces the 


APPENDIX 

A. SYSTEMATIC UNCERTAINTIES IN SFR TRACERS 

The effect of a diffuse component in the Xsfr measure¬ 
ments used in our grid analysis is quantified by using a simple 
approach. Our approach assumes that the diffuse component 
not associated to SF is mainly dominated by a cirrus com¬ 
ponent in the 24/im map. Additionally, our approach assumes 
that the diffuse emission is negligible in regions of high SF, so 
the Xsfr measurements from 24/im+FUV is a good estimate 
for the true Xsfr in this regime. On the other hand, it assumes 
that the Xsfr measurements from Ha are less affected by a 
cirrus component at low SF, so any deviation from a linear re¬ 
lation between Xsfr from 24/zm+FUV and Xsfr from Ha is 
assumed to be due to a cirrus component in the 24/im map. 

Our approach is as follows. First, we smooth the 24//m map 
using a gaussian kernel of a given size. Then, the smoothed 
map is multiplied by a factor that accounts for the fraction 
of the emission associated to the diffuse component. Finally, 
this scaled and smoothed map is subtracted from the original 
24//m map. Thus, our approach relies on two parameters: the 
size of the gaussian kernel, 0ker, and the assumed fraction of 
diffuse emission in the smoothed 24/jm map, / DE . We have 


TABLE 8 

Pearson correlation coefficient for Ssfr.hq vs. SsFR, 24 /jm+Fuv 

LINEAR CORRELATIONS. 


0ker 


No diffuse 



/de 





correction 

0.2 

0.3 

0.4 

0.5 

0.6 




NGC 

6946 




225" 

rp 

0.912 

0.926 

0.929 

0.920 

0.864 

0.907 

NGC 628 

170" 

rp 

0.849 

0.862 

0.869 

0.876 

0.882 

0.888 

M101 

170" 

rp 

0.928 

0.932 

0.930 

0.927 

0.916 

0.885 


highest value of r P =0.929. For NGC 628, the highest value of 
r P =0.888 is obtained with / DE =0.6. Finally, / DE =0.2 gives the 
highest value of r P =0.932 for M101. 

In Figure [25] we illustrate the effect of applying the 24//m 

cirrus correction to the Xsfr measurements for NGC 6946. In 

—1 —2 

the left panel, we observe that at Xsfr ~ 0.01 M 0 yr kpc ", 
the measurements of Xsfr from 24/mi+FUV overestimates 
the Xsfr from Ha. After a cirrus component of 24//in is 
subtracted using /de=0.3 and d^ er =225'', both approaches 
to estimate Xsfr present a tighter correlation for Xsfr < 
0.01 Mq yr" 1 kpc -2 . 
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Fig. 25.— Left: Comparison between Ejfr using Ho and the Sjfr using FUV + 24//m for NGC 6946 (same as Figure[7]but for the uniform grid). Black 
dots represent on-arm regions, while blue dots correspond to inter-arm regions. Right: Same as left panel but now a subtraction of a 24//m smooth component 
is applied to Esfr using FUV + 24//m. As in FigurelT] we use Ah q = 1.0 magnitude extinction for regions in NGC 6946. After the 24//m smooth component 
subtraction is applied, the linear relation between the two Dsfr approaches is better recovered. 


Karachentsev, I. D., Karachentseva, V. E.. Huchtmeier, W. K., & Makarov, 
D.I. 2004, AJ, 127, 2031 
Kelly, B. C. 2007, ApJ, 665, 1489 
Kennicutt, R. C„ & Evans, N. J. 2012, ARA&A, 50, 531 
Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 

Kennicutt, Jr., R. C., Armus, L., Bendo, G., et al. 2003, PASP, 115, 928 
Kennicutt, Jr., R. C., Calzetti, D., Walter. F., et al. 2007, ApJ, 671, 333 
Kim, W.-T., & Ostriker, E. C. 2002, ApJ, 570, 132 
Koda, J„ Sawada, T„ Wright, M. C. H„ et al. 2011, ApJS, 193, 19 
Koda, J., Scoville, N.. Hasegawa, T., et al. 2012, ApJ, 761, 41 
Kravtsov, A. V. 2003, ApJ, 590, LI 

Kruijssen, J. M. D„ & Longmore, S. N. 2014, MNRAS, 439, 3239 
Kuno, N„ Sato, N„ Nakanishi, H„ et al. 2007, PASJ, 59, 117 
Lada, C. J., Forbrich, J., Lombardi, M., & Alves, J. F. 2012, ApJ, 745. 190 
Larson, R. B. 1981, MNRAS, 194, 809 
Leroy, A. K„ Walter, F„ Brinks, E„ et al. 2008, AJ, 136, 2782 
Leroy, A. K„ Walter, F„ Bigiel, F„ et al. 2009, AJ, 137, 4670 
Leroy, A. K., Bolatto, A., Gordon, K., et al. 2011. ApJ, 737, 12 
Leroy, A. K„ Bigiel, F„ de Blok, W. J. G., et al. 2012, AJ, 144, 3 
Leroy, A. K., Walter, F., Sandstrom, K., et al. 2013, AJ, 146, 19 
Liu, G., Koda, J., Calzetti, D., Fukuhara, M., & Momose, R. 2011, ApJ, 735, 
63 

Louie, M„ Koda, J., & Egusa, F. 2013, ApJ, 763, 94 
Mandel, K. S., Narayan, G„ & Kirshner, R. P. 2011, ApJ, 731, 120 
Markwardt, C. B. 2009, in Astronomical Society of the Pacific Conference 
Series, Vol. 411, Astronomical Data Analysis Software and Systems 
XVIII. ed. D. A. Bohlender, D. Durand, & P. Dowler, 251 
McKee, C. F„ & Ostriker, E. C. 2007, ARA&A, 45, 565 
Momose, R., Okumura, S. K., Koda, J., & Sawada, T. 2010, ApJ, 721, 383 
Momose, R., Koda, J., Kennicutt, Jr., R. C., et al. 2013, ApJ, 772, L13 
Narayanan, D., Krumholz, M. R.. Ostriker, E. C., & Hernquist, L. 2012, 
MNRAS, 421,3127 


Ossenkopf, V., & Mac Low, M.-M. 2002, A&A, 390, 307 
Passot, T., & Vazquez-Semadeni, E. 1998, Phys. Rev. E, 58, 4501 
Rahman, N., Bolatto, A. D., Wong, T., et al. 2011, ApJ, 730, 72 
Rebolledo, D., Wong, T., Leroy, A., Koda, J., & Donovan Meyer, J. 2012, 
ApJ, 757, 155 

Rosolowsky, E., & Leroy, A. 2006, PASP. 118, 590 
Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999, ApJS, 
124, 403 

Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5 
Schinnerer, E., Meidt, S. E., Pety, J., et al. 2013, ApJ, 779, 42 
Schlegel, D. J.. Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schruba, A., Leroy, A. K., Walter, F., Sandstrom, K., & Rosolowsky, E. 
2010, ApJ, 722, 1699 

Schruba, A., Leroy, A. K„ Walter, F„ et al. 2012, AJ, 143, 138 
Sharina, M. E., Karachentsev, I. D„ & Tikhonov, N. A. 1996, A&AS, 119, 
499 

Sheth, K„ Vogel, S. N„ Wilson, C. D„ & Dame, T. M. 2008, ApJ, 675, 330 
Shetty, R., Beaumont, C. N.. Burton, M. G., Kelly, B. C., & Klessen, R. S. 
2012, MNRAS, 425, 720 

Shetty, R„ Kelly, B. C., & Bigiel. F. 2013, MNRAS, 430, 288 

Solomon, P. M„ Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 

Tully, R. B. 1988, Journal of the British Astronomical Association, 98, 316 

Wada, K„ Baba, J., & Saitoh, T. R. 2011, ApJ, 735, 1 

Wada, K„ & Koda, J. 2004, MNRAS, 349, 270 

Wada, K„ & Norman, C. A. 2001, ApJ, 547, 172 

Walter, F„ Brinks, E„ de Blok, W. J. G„ et al. 2008, AJ, 136, 2563 

Wong, T„ Hughes, A., Ott, J., et al. 2011, ApJS, 197, 16 

Wyder, T. K.. Martin, D. C., Schiminovich, D., et al. 2007, ApJS, 173, 293 



























